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FOREWORD 


This report (Volume I) is a technical summary of the effort on Modifica- 
tion 3 of Contract NAS 8-27228 carried out between May 18, 1972 and November 
21, 1972. The contract was begun on August 25, 1971 as a continuation of the 
technological development efforts carried out under Contracts No. NAS 8-24683 
and NAS 8-26157. On September 14, 1971 the work was verbally redirected by 
Marshall Space Flight Center to cover definition, design, manufacture, test 
and flight qualification of an experiment package which would permit the 
carrying out of zero gravity supercooling and nucleation experiments in the 
Skylab mission as an extension of the M512-related experiments. The objectives 
of the supercooling and nucleation experiments are described in an Experiment 
Proposal dated June 18, 1971 submitted to Marshall Space Flight Center by 
Dr. R. T. Frost and Professor D. J. Turnbull as Principal Investigators. 

This document served as the initial guide for the redirected effort. 

On December 10, 1971 Marshall Space Flight Center approved a number 
of specific tasks (Modification 1) which served as a more detailed framework 
for the new redirection. Work was greatly accelerated in order to meet the 
very stringent deadlines set by the M512 qualification and launch schedule. 

In February 1972 NASA made a decision to terminate the accelerated 
effort aimed at furnishing a supercooling and nucleation experiment for the 
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Skylab and issued a stop work order on February 18. On May 18, 1972, a 
completion effort was defined and mutually agreed upon (Modification 3) which 
would result in preserving those parts of the technical effort which would be 
most useful with regard to further technology developments in the zero gravity 
melting and solidification of metals and semiconductors. Since a large amount 
of the technical effort had been determined by the specific compatibility require 
ments with the existing M512 experiment apparatus and the Skylab interfaces, 
this termination effort was primarily directed to documenting and preserving 
technology and hardware which would be of most general interest for future 
experiment developments rather than documentation of the details of the 
technical effort most specifically related to the M512 and Skylab constraints. 

To this end, attention was concentrated upon assembly and delivery of labora- 
tory breadboard hardware consisting of a working device for automatic position 
sensing and control of freely floating diamagnetic spheres compatible with the 
Marshall Space Flight Center 400 foot drop tower test program. The new 
documentation of the various phases of the engineering work was limited to 
only those few items of most general interest for further work in this new 
technology area. 

In Volume I, brief summaries are given of the selected technical areas 
as well as descriptions of the breadboard hardware and drop test package which 
is being delivered as part of this effort. Volume II of the report contains the 
documentation on the Modification 1 contract work specifically related to the 
Skylab experiment as it existed at the time of contract completion. 
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1.0 INTRODUCTION 


The material which follows describes briefly the technology development 
done in conjunction with defining a facility to be used on the Skylab mission for 
electromagnetic suspension of small, molten spheres in the weightless space 
environment. This facility was planned as an add-on for the M553 sphere - 
forming experiment. That experiment utilizes an electron beam for melting 
the specimen and it was planned to use the same source of energy for the 
electromagnetic suspension demonstration. Electromagnetic suspension and 
position was to be employed to control the position of the specimens so that 
they would not contact any other materials or objects during solidification. 

In the proposed add-on system, electromagnetic forces are developed 
by the interaction of the fields from a system of current-carrying coils, with 
the dipole moment induced by this same coil system in the floating melted 
specimen. Position and velocity of the levitated spherule are continuously 
monitored by an optical position sensor. A temperature history of the melt 
during cooling and solidification was to be provided by means of a simple 
radiation pyrometer measuring the ratio of radiance in two infrared wave 
bands. A position control system introduces position and position rate infor- 
mation into a servo loop to control the current levels in several nearby mag- 
netic coils. Although positioning forces as high as 50 to 100 dynes would be 
utilized for initial deployment of the spherule at the instant of melting and for 
recovery of the specimen after solidification, it was expected that position 
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and velocity errors would be nulled within the first few seconds after melting 

so that only miniscule forces would have to be applied to this spheroid during 

solidification, provided that solidification did not begin for several seconds 

after melting. "Miniscule forces" here refers to those forces required to 

overcome residual accelerations of the Skylab facility which are expected to 
-5 

be in the 10 g range. 

1.1 GENERAL DESCRIPTION 

The Electromagnetic Levitation System (EMLS) was to utilize the 
spherical work chamber of the sphere -forming experiment. Figure 1.1-1 
illustrates, schematically, the various elements of the proposed system posi- 
tioned in this chamber and includes a perspective view of the specimen wheel 
with the electromagnetic positioning coils, position sensors and control elec- 
tronics. The electron beam of the facility enters the sphere at the lower left 
end as indicated by a broken line. Also shown by a broken line is a small 
cylinder which was to serve to capture the specimens after solidification. A 
detailed description of the Electromagnetic Levitation System as planned for 
the Skylab is given in the Experiment Implementation Plan, Section A of 
Volume II of this report and will not be repeated here. 

1.2 BACKGROUND 

Under Contract NAS 8-24683 analyses and experiments were carried out 
on a position control system consisting of six coils mounted as three orthogonal 
pairs. The coils consisted of typically hundreds of turns of relatively small 
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Figure 1. 1-1. EMLS in M512 Vacuum Chamber 
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diameter wire and were operated at frequencies in the 10 kHz region. The 
coils could be considered to lie on the faces of a cube circumscribing the 
specimen position. The position of the conducting sphere relative to the center 
of the cube was monitored by means of the change in inductance of the individual 
coils due to the reflected impedance of the eddy current paths in the spherical 
specimen. The position error signals, consisting of the relative changes in 
impedance of the six coils, was fed into a servo system which automatically 
regulated the relative power inputs to the six coils in such a manner as to urge 
the spherical specimen to the center of the cube and to damp oscillatory motions. 

Under Contract NAS 8-26157, this concept was developed into a zero 
gravity position and velocity control servo which was tested in the MSFC 400 
foot drop test facility and in the KC135 USAF ballistic flight program. This 
package incorporated electro-optical sensors for determining the position 
errors and error rates for the freely floating specimen. Six coils were again 
used for position control and simultaneously served as the inductor in resonant 
tuned circuits to improve power coupling efficiency at operating frequencies in 
the 100 kHz regime. The electro-optical position sensing development was 
undertaken to widen the range of specimen resistivities which could be handled 
by the facility since the electromagnetic position sensing scheme is restricted 
to relatively good conductors. The use of the coils simultaneously as producers 
of the magnetic field for position control and as tank circuit inductors resulted 
in only relatively weak available electromagnetic forces because of the high 
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tank inductance required for use with available capacitors, causing the coil 
configuration to be solenoids having appreciable length. 

1.3 SUMMARY 

At the beginning of the contract work reported here, it was decided to 
separate the field producing and tank circuit inductor functions by use of ferrite 
core transformers which coupled a resonant tank circuit to the field producing 
coils. It was also decided to explore the relative merits of employing only 
four position control coils arranged so as to lie on the faces of an imaginary 
tetrahedron. This configuration has the advantage of giving a somewhat more 
open structure to the coil assembly and reduce the number of power amplifiers 
required. Mock-ups and tests as well as analyses were carried out and a 
comparison with the six coil configuration was begun. When the contract was 
redirected towards furnishing a flight model equipment for the Skylab Super- 
cooling and Nucleation experiment it was decided that the more open structure 
of the tetrahedron as well as a saving in number of power amplifiers was more 
compatible with the constraints of the already existing M553 specimen changing 
wheel, electron beam heater and vacuum chamber. For that reason most of 
the work directed towards the Skylab experiment application was based upon 
the assumption of a four coil tetrahedral configuration. The tetrahedral con- 
figuration required implementation of coordinate converter circuitry since 
there exist only three independent coordinates describing the position error of 
the specimen, whereas excitation of the separate coils results in force components 
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lying in the four non-orthogonal directions described by the coil axes. The 
requirement for this additional signal circuitry was considered to be a dis- 
advantage outweighed by the superiority in terms of coil geometry and number 
of high powered circuits and coil leads as compared to the six coil system. 

At the time the Skylab Supercooling and Nucleation experiment effort was 
terminated, the coordinate conversion signal processing electronics had not 
passed the conceptual design stage. In the interest of furnishing a working 
drop tower demonstration package for future experiment development in this 
area on the limited funds available, it was determined that a six coil unit could 
be assembled more readily because of the relatively large expense that would 
be entailed in detailed design and construction of the coordinate converter 
electronics. It was also considered that the advantages of coil configuration 
and reduced number of power amplifiers were associated primarily with the 
constraints determined by the existing Skylab equipment and would probably 
not outweigh the disadvantage of the coordinate converter electronics for future 
experiment development unconstrained by details such as the M553 specimen 
selector wheel. A comparison of the four-coil and six-coil configurations is 
considered in Section 2.2 of this report. 

Because of the necessity for performing certain basic tests on the electro 
magnetic force field configuration with the four-coil system on the accelerated 
program aimed at the Skylab opportunity, it was necessary to carry out tests 
on a four-coil servo mock-up at a time before the coordinate converter 



electronics was to be available. To this end, a demonstration servo controlling 
the position of a neutrally buoyant sphere floating in liquid of controlled density 
was carried out with the four-coil configuration. A simplified position error 
sensing system and position control servo without coordinate converter was 
assembled for these tests. Although the position sensing and control servo 
circuits were somewhat different than those adopted as the reference design for 
the Skylab experiment, it was agreed to deliver this working four-coil position 
servo as part of the Mod 3 termination effort in order to preserve the four -coil 
breadboard hardware development which had been carried out to that date. 

The position error sensor which was adopted as the reference design for 
the Skylab Supercooling and Nucleation experiment consisted of two separate 
electro-optical sensors viewing the specimen in two orthogonal directions. 

Each sensor incorporated a quadrant silicon detector which could determine 
the specimen position transverse to the viewing axis. With the two detectors 
the three orthogonal position errors of the specimen could be measured with 
some redundancy. It was decided to incorporate these detectors into the 
deliverable position control drop tower demonstration package as the most 
valuable way of preserving this important technology development. These 
detectors are described in Section 2. 6 and are also referred to in the drop 
tower package description in Section 3.0. 

For monitoring of the specimen temperature during solidification, a 
miniature two-color pyrometer was designed and breadboarded. The decision 
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to utilize a two-color optical pyrometer rather than a total radiation pyrometer 
was dictated by the desire to observe specimen temperature changes in the 
absence of detailed knowledge of specimen emissivity, particularly during 
phase transformations at which discontinuities might occur in emissivity. In 
the Skylab Supercooling and Nucleation experiment it was planned to incorporate 
the pyrometer into the same housing as the optical position sensors. A working 
breadboard miniature two-color pyrometer is being delivered under the terms 
of Mod 3 of this contract but was not incorporated into the drop tower package, 
which at present is designed only for position control demonstrations with cold 
solid specimens. 

Because of the absence of laboratory test data on vaporization rates from 
the molten specimens in the vacuum provided by the existing Skylab facilities 
(particularly because some experiment specimens were not determined at the 
time), there was some concern for problems arising from plating of the optical 
position sensors by metal vapors. This effect was to be minimized by an 
arrangement by which the solid state optical sensor element viewed the hot 
specimen in reflection by a front service mirror. Baffles were considered 
which would shield the detector from vaporization products traveling directly 
from the molten specimen in vacuum. Because of the possibility of detector 
contamination by specimens which might be proposed by other investigators to 
which the detector would be especially susceptible (it was considered that the 
detector could survive a wide range of specimen experiments), a back-up 
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position control system was considered utilizing the electromagnetic detection 
principle originally studied under Contract NAS 8-26157. This work is dis- 
cussed briefly in Section 2.5. 1. 

An analysis is given of the possibilities for specimen rotation control 
using the electromagnetic position control fields. This is important with respect 
to inhibiting or controlling rotation rates induced by single -phase induction 
motor action in freely floating conducting liquid masses. In addition, possi- 
bilities exist for shape forming masses into oblate spheroids by inducing 
controlled rotation rates by means of a two -phase induction motor action which 
can be arranged by proper phasing of the position control coils. 

Section 3.0 gives some detailed description of the demonstration drop 
tower package for electromagnetic position control. This package was designed 
to be compatible with the MSFC 400 foot drop tower facility. Diagrams are 
given showing connections for power and telemetry which should be fully com- 
patible with the carriage which was developed at MSFC for use in both the drop 
tower and the KC135 aircraft. 

Listings of the computer programs which were written for the magnetic 
field and force analyses are given in an appendix. 

In the course of performing the work described in this report, numerous 
drawings were generated. Some of these drawings are included in this document 
but many were not of a nature to be referenced here. However, a list of all 
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drawings prepared is included in Volume II. Microfilm copies of all drawings 
have been provided to the Contracting Officer 1 s Representative and prints may 
be obtained from him, if desired. 
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2. 0 TECHNOLOGY DEVELOPMENTS 


In this section the technology developments are described which are of 
most general interest for applications to subsequent efforts in the new area of 
containerless processing of metals and semiconductors. 

2. 1 FOUR COIL OPTIMIZATION 

In this section we will discuss the fourrcoil tetrahedral arrangement and 
its optimization. 

The function of the four coils in the tetrahedral configuration is to provide 
the force, within the available power, required for specimen deployment or 
retrieval, damping, positioning and, after solidification, ejection or stowage. 
This force is obtained by selective excitation of one or more coils by radio 
frequency currents controlled by the position sensing subsystem. The theory 
describing this force is treated in Section 2 of the final report dated June 15, 
1971 on Contract NAS 8-26157. It was shown there that force is dependent on 
the following parameters. 

1. current magnitude in driven coil 

2. driven coil size and shape 

3. specimen size 

4. specimen to coil distance 

5. skin depth of specimen. 
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Efforts were made to optimize the force capability of the coils under the 
following constraints, 

A. Mechanical 

• M512 vacuum chamber size 

• Specimen size 

• Electron Beam Clearance 

• Specimen, Specimen Wheel and Indexing Motor Clearance 

• Specimen Stowage Path. 

B„ Optical 

• Position Sensor Field of View 

• Temperature Sensor Field of View 

• Camera Field of View 

• Visual Field of View 

C. Electrical 

• Input R. F. Power (Thermal) 

• Coil Interactions 

• High Voltage 

A four-coil facility was identified that was compatible with these con- 
straints and which could handle specimens up to 1 cm in diameter. The coils 
were 2 cm inside diameter, spaced tangent to an inscribed sphere of 4 cm 
diameter. These coils were each wound with 3 turns and had a cross section 
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approximately 0. 4 cm. The coil and its leads are cut from one piece of copper 
strip about 0. 08 cm thick and is self-supporting. 

Referring to the expression* for the total force exerted on the sphere by 
the applied field (Equation 2-4 of the referenced final report): 

3 dH 

F = - 277 R^ /LiH — G(X) where 
= Sphere Radius, 

G(X) = Body Force Function, discussed in more detail below, 

Z = Axial Distance from Coil to Sphere, 
we can describe how each parameter can be varied to optimize force. 

2. 1. 1 Sphere Radius (R^) 

The force of a given field and field gradient varies as the cube of R^. For 
maximum force we see that R^ must be kept as large as possible. We will see 
later where G(X) also requires a large R 2 for appreciable forces. For free 


* In the formulae in this section, the MKS system of units is used, i. e. meters 
kilograms and seconds. The field intensity is expressed in Maxwells. In 
some later sections cgs units are used, i.e. centimeters, grams and seconds, 
because of their greater convenience for apparatus of the dimensions con- 
sidered here and because of the convention of expressing surface tensions and 
viscosities in these units. Because of the non-triviality of conversion 
between these sets of units when one is considering magnetic field strengths 
and induction, a short discussion of units and their inte r -conversions is given 
in Appendix C. We will encounter mixed units and the necessity for conversion 
in Section 2.7. 
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body acceleration forces (g forces) the mass varies as the cube of R so R is 

2 2 

not significant for acceleration except, as noted above, in G(X). In the event 


there is wetting of the pedestal holding the molten specimen it is important to 
provide adequate electromagnetic faces to overcome surface tension forces 


binding the specimen to the sting on which is is melted, as was the arrangement 
using the M553 specimen wheel. The sphere radius R then must be kept as 
large as practical and the coil system must have a free volume sufficient to 


handle a sphere of this size. 



where r is the slant distance to the winding and i is the current in amperes. 

The force exerted on the specimen is proportional to H and therefore varies 
as the square of the coil radius, the inverse cube of the coil-specimen distance, 
and directly with current i. Immediately we have contradictory requirement 
because increasing a makes the coil assembly larger and that makes r larger 
for a specimen located at the center of the coil assembly. Because H is more 
sensitive to r than to a, some advantage is gained by keeping the coil assembly 
compact. 
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2. 1. 3 Magnetic Field Gradient 

In the center of the coil the field H is greatest but the change in field, 
dH 

-j-g- , is zero. At a far distance along the axis, the field and field gradient 
are small. As we move away from the coil the field starts to reduce and the 
gradient increases rapidly. The force, which is proportional to the product of 
the field and the field gradient, must therefore increase as we proceed away 
from the coil plane, but vanishes at large distances. It is thus obvious a 
maximum force exists somewhere. This maximum is found to occur at a 
distance from the coil of approximately one -third the coil diameter. 

2.1.4 Body Force Function, G(X) 

R 2 

Figure 2.2-1 is a plot of the function G(X). X = -g— where is the 
specimen radius as before and 6 is the conducting skin depth. Nothing can 
be done to the coil to optimize G(X). The specimen size can be maximized 
as mentioned previously but will affect the volume of the coil assembly. 

Because X needs to be at least 3, so that G(X) -* 1, the skin depth must be 
kept shallow. Therefore, for any specimen size and resistivity, there is 
a minimum frequency which must be exceeded in order to yield a G(X) near 
unity. 

2.1.5 Selection of an Operating Frequency 

The frequency of operation is significant because it must be above that 
required to make G(X) ~ 1 for the highest resistivity material being considered 
for the specimen. For R^ = 0. 5 cm, 6 ^ 0. 17 cm, and a specimen material 



2-5 







with a resistivity of 100 fiQ -cm, this requirement dictates a minimum fre- 
quency of 900 kHz. 

We have described how parametric variation affects force and now must 

consider the constraint of available power. Changing parameters to maximize 

force alone is insufficient; maximum force per unit power input is the desired 

criteria. The major R. F. power loss of the system occurs in the coils. In 

our consideration of the magnetic field we saw that the field is proportional 

to the current, the coil radius squared and the inverse cube of the specimen 

distance. The coil loss is proportional to the resistance and therefore to the 

2 

coil radius, a, and the field is proportional to a . We gain in force per unit 
power by using the largest a practical. The need for clearances, however, 
puts a very definite limit on the coil radius, a. The coil loss is proportional 
to coil resistance so we find that the higher the frequency the greater the loss 
due to the decrease in conducting skin depth of the coil material. As a minimum, 
the frequency should not go above that which produces a skin depth of half the 
coil conductor thickness. For copper at 90 kHz, the skin depth, 6, is 0. 022 cm. 
With a coil conductor thickness of 0. 08 cm, 6 should be in the order of 0. 04 cm 
or larger. This condition would occur at 27 kHz. Above this frequency, the 
coil loss increases by the square root of the frequency. However, we have 
found that a minimum operating frequency approaching 90 kHz must be used 
because G(X) drops off below much faster rate than the square root of the fre- 
quency below 90 kHz. Therefore the optimum frequency for a 100 /i r-cm 
specimen is in the neighborhood of 90 kHz. 
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From this discussion, we can conclude that the maximum force per unit 
power is always lower for the higher resistivity materials in any given coil 
configuration. 

The previous analysis indicates that the copper conductors can be thinner 
and the number of turns increased. Such a coil is feasible for this application 
if the temperature rise in the coils for full power input is within the limitations 
of any insulators needed to survive grazing contacts with the hot specimen. 

2.1.6 Impedance Matching Transformer 

Another item which entered into the optimization of the coils was the 
impedance matching transformer. This unit had a one turn secondary which 
for the available core being used could drive a coil of at most three turns. For 
a four turn positioning coil, a second secondary turn was needed at the power 
level used. This would double the turns in primary winding. Overall, the 
additional loss in the impedance matching transformer with the larger number 
of turns would offset the gain in going to four turns in the positioning coil. 

A four turn coil is feasible if a special (non-standard size) core is specified 
for the impedance matching transformer, a procedure that could be used on 
any flight program. 

2.1.7 Mechanical Configuration 

In the simplest geometric arrangement of the four-coil system, each 
coil would be oriented to be in the face of a simple tetrahedron as shown in 
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Figure 2. 1-2. This initial arrangement permitted the specimens to enter the 
coil system along path A, the electron beam to enter via path B (which also 
includes a continuing clear exit path in case the beam is ever permitted to 
progress that far) and the ejected specimens to leave via C (or in the direction 
opposite to that of C for "on-the -wheel" stowage). This initial concept was 
considered to be unsatisfactory primarily for the following two reasons. 


A. If a specimen should stick to its support on the wheel after having 
been melted, the forces which may be created at the location of the 
support may be too small to push the specimen off the support. 

B. If the specimen should readily free itself from the support and push 
itself away from it - as a consequence of it assuming a spherical 
shape and of non- wetting of the support - the forces on the opposite 
side of the coil system, midway between two coils may be insufficient 
to prevent the specimen from leaving the coil system prematurely. 
Analysis of this situation is included in Section 2. 3 below. 


It was essential to the success of the experiment that the specimen come 
free of its support and then remain within the coil system without touching the 
coils until it again becomes solid. Since no useful data on the release of typical 
molten specimens from supports in a freely falling experiment were available 
during design of the coil system, it was decided that the set of coils, and 
possibly the specimen with its support, must be rotated to place the specimen 
at a location at which the forces were known to be larger. This resulted in a 
tetrahedronal coil configuration and wheel modification shown in Figure 2. 1-3. 
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Figure 2. 1-3 Modified Tetrahedronal Coil Arrangement 




2. 2 FOUR-COIL VERSUS SIX-COIL CONFIGURATION COMPARISON 
2. 2 0 1 Size 

For three dimensional control the four-coil configuration requires only 
two thirds the number of high power components of a comparable six-coil con- 
figuration. Because the power components are the largest parts used, it 
follows that the four coil system had the best chance to fit within the small 
space available for this experiment. 

The small signal parts count however will rise significantly. These parts 
are used to perform the coordinate conversion and the servo force vector 
resolution. This results from the four-fold symmetry of the tetrahedral con- 
figuration as compared to the existence of only three independent coordinates 
for the position errors. The position sensors are on three corresponding 
orthogonal axes. In general, output from three coils is needed to generate force 
vectors that are parallel to these axes. A further discussion of the tetragonal 
coordinate system, coordinate converter and the servo is found in Section 2. 3. 

In the six-coil system, no coordinate converter is required. In this con- 
figuration, axes of the coils define three orthogonal position measurement axes, 
and the force vectors can easily be derived as proportional to the negative 
of the position error vectors. Each of the three sensing systems controls only 
the two coils along that coordinate axis and has no coupling affect into the other 
coil systems. In total piece part count the six-coil system has fewer parts than 
the four-coil system. The size and weight of the six-coil system is greater by 
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about twenty percent than those of the four-coil system because of the larger 
count of high power components. 

2. 2. 2 Force Considerations 

When the force vector created is along a coil axis and directed away from 

that coil, relatively high force levels are achieved. When the force is directed 
at a significant angle off-axis, we find the force greatly reduced. In fact, if 

the specimen is near a coil and on-axis (in a system as defined in Section 
2. 3. 2. 3), it could experience over 50 dynes of force. On the other hand, if the 
specimen is approaching a corner of the four coil system (the part of the tetrahedron 
on a coil axis but beyond the center of the coil assembly) the force oo uld not 
exceed about 5 dynes at the very modest total power available from existing M512 
circuitry. Further, if the specimen moves 1 cm or more from the center toward 
a corner of the four-coil system, the force is zero or negative and the specimen 
can be lost by being ejected from the enclosure. The difficulty presented to the 
servo by this gross non-symmetry of available forces in the four coil system was 
being resolted at the termination of work on this contract as reported in Section 
2. 3. 3 below. In view of the forces resulting from the specimen either wetting 
the substrate or being repelled by it as the molten specimen assumes a spherical 

shape, and the lack of experimental data in this area, considerable uncertainties 
remained. 

The forces from a given coil in the six-coil system are no greater than 
those from an identical coil in a four-coil system. The six-coil system, how- 
ever, can contain a specimen so long as the specimen is within 2 cm of the 
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enclosure center. Further, if the specimen were near a coil and directed 
toward a corner (where the force in a six-coil system is no greater than in a 
four-coil system), it must pass near two coils both of which will direct it 
toward the center thereby preventing the specimen from approaching the 
corner and reducing the probability of it being ejected out of the enclosure. 

A further operational advantage of the six-coil system is that for any 
force vector directed to the center region of the enclosure, an equal and 
opposite vector can be generated after the specimen passes the center. Where 
there was gross non- symmetry in the four-coil system, we find symmetry in 
the six-coil system and the servo difficulty noted previously in the four-coil 
system does not exist. 

2. 2. 3 Interface Effects 

Because the EMLS wheel was to have been attached to the M553 Sphere 
Forming Experiment' s motor shaft and because the distance between this shaft 
and the electron beam had been determined, the diameter of the specimen wheel 
was required to remain unchanged from that for the M553 experiment. The 
specimen size was selected to be the same as that used in M553, that is, a 
one centimeter diameter sphere when molten, with consideration for some to 
be of a smaller size to accommodate metallurgists whose experimental tech- 
niques require the use of small quantities of material. The use of a smaller 
specimen in a system designed to process larger specimens poses a possible 
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problem to the position control system because as the ratio of sphere size to 
coil size decreases, the forces which can be exerted upon the sphere decrease 
as discussed earlier in Section 2. 1. 1. 

Because of concern of the dynamics of a specimen within the coil system 
and the possibility that a specimen might be lost through a weak corner (i.e. , 
a region in which the force the coils can exert upon a specimen is relatively 
small) as discussed in Section 2. 3. 3, a six-coil system with nearly square 
coils was considered as shown in Figure 2. 2-1. Note that one of these coils 
had to be distorted somewhat to provide a clear entry for the incoming 
specimens. For comparable volumes of specimen confinement within the 
coils the two coil systems were found to produce forces of comparable sizes. 

It had been found from previous contract work that the six-coil system gave 
suitable position control and damping from experimental tests and analyses. 

No work was done, however, under the present contract to analyze specimen 
Wetting to the sting except for initial conceptual theory. 
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Figure 2.2-1. Possible Six Coil Configuration 


2. 3 FOUR-COIL POSITION CONTROL SERVO 


2. 3„ 1 Introduction 

The four-coil tetrahedral configuration was examined with respect to 
specimen control. Mathematical equations defining the ideal control process 
were derived. These equations were then integrated with and modified by the 
results of theoretical and laboratory performance analysis of the actual coils. 

On one hand, the defining control equations had to be kept as simple as 
possible in anticipation of their ultimate implementation via the servo elec- 
tronics. On the other hand, the actual field forces produced by one coil or by 
two or three coils acting simultaneously are quite complex and are discussed 
in Section 2. 3. 2.4. The exact solution of the system of resulting differential 
equations is however formidable and implementation of the solution in a servo 
control requires approximations. 

The approach taken was minimization of the servo electronics by 
simplifying assumptions while maintaining control over a specimen contained 
in the complex coil force field. This objective was met. The details of the 
computations and simulations that have been performed are presented below. 

2. 3. 2 Derivation of Servo Control Equations 

2. 3. 2. 1 Field Approximation and Comparison to Measurements 

The sequence of events as they appeared in the servo control operation 
is as follows: 
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1. Specimen position is sensed and velocity calculated by the position 
sensors. 

2o A coordinate transformation is made to the defined axes. 

3. Based on the spring-mass -dashpot equations of motion, a vector 
force is calculated which, if applied to the specimens, would drive 
it toward the tetrahedron origin (0, 0, 0). 

4. The calculated vector force of step 3, F , is resolved into vector com- 
ponents for each of the four coils. 

5. An equation is solved, applicable to each coil, which relates the force 
desired from that coil and the input current required to produce that 
force. 

6. The current of step 5, is delivered to the coils ^.nd the coils respond by 
exerting a total vector force on the specimen, F^,. 

7. After some At, say, 0.01 sec., 1. begins again. As the specimen 
nears (0, 0 # 0) damping is accomplished and the specimen eventually 
comes to rest at (0, 0, 0). 


In the highly idealized system the command force, F^, is equal to the 
delivered force, FJJ,. If the model assumed to compute the equations of Step 3. 
posses a stable solution, then with F^, = F,^,, the specimen will be precisely 
controlled and be brought directly to (0, 0, 0). i n practice, however, taking 
into account all the field interactions would require a very complicated control 
system. 

With some reflection, it is clear that the equation of motion assumed in 
step 3. above is somewhat arbitrary. All that is required is that these equations 
generate commands for F,^ appropriate to the desired system operation; strong 
forces on the specimen, directed toward the origin, when the specimen is far 
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(~ 1 cm) from the origin; forces directed opposite to velocity and proportional 
to velocity (damping); and zero force when the specimen is at (0, 0, 0) and at 
rest with respect to the control volume. 

A tractable equation of motion is provided by a spring-mass -dash pot 
analogue. The solution of these equations can be implemented by the four-coil 
system for suitably chosen values of the spring constant, k. The controlling 
equations are of the form 

F_ = m = -kp + b p- (1) 

T , 2 ^ dt 


Where: F^, is the force applied to the specimen and whose mathematically 

derived direction is toward the origin of the tetrahedral system (0, 0, 0); p is 

the displacement vector of the specimen from (0, 0, 0); k and b are the ' spring 

1/2 

constant' and damping term, respectively; b = (4 m k) for critical damping; 
and m the mass of the specimen. 


For a specified value of k, the controlling equations are defined for 
critical damping. The value of k selected depends primarily on the capability 
of the individual coils, the resultant field and the specimen material. 

The value of k must be judicially chosen. Too high a value would result 
in the coils operating near or at maximum current until the specimen gets very 
near (0, 0, 0). Too low a value of k would result in the neglect of the system' s 
full capability. Ideally, a unique value for k would be use with a given coil 
geometry for each specimen. 
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2. 3, 2. 2 Current to the Coils 


Equation (1) specifies the command force F , i.e., the force to be applied 
to the specimen. Whether or not the specimen actually experiences the F 

depends upon the equation used to relate the command force to the current to 
be delivered to the coils and even more so upon the actual field behavior of the 
forces generated. 

The solution of Equation (1) assumes that the forces generated off-axis of 
the coil are identically equal to the force that would be calculated, for the same 
distance from the coil, on the center line. That is, Equation (1) contains no 
provisions for variation of the field over the radius of the coil. There is no 
simple procedure by which the off-axis effects can be handled mathematically, 
save by a solution of the equation of motion with these effects included. This 
would lead, however, to a complex servo electronics and, as will be shown, 
such complications are necessary. 

The block diagram of Figure 2. 3-1 outlines the ideal operation of the 
system and the approximate approach for which the electronic design has been 
completed. 

As noted on Figure 2.3-1, the ideal servo analysis would include the 
electromagnetic field and force expressions derived by Fromm and Jehn*. 


* Fromm, E. and Jehn, H. , "Electromagnetic Forces and Power Absorption 
in Levitation Melting," Brit. J. Appl. Phys. , 1965, 16, 653. 
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Here the command forces are equal to those delivered, F^, = F^ , where the 
i subscript denotes the acting coil. 

Equation (1) yields the desired electronically implemented command forces. 

Employing the approximations to the electromagnetic force equations at 4b of 

Figure 2.3-1 yields currents that, by definition, yield the command forces. 

However, if the specimen is off-axis of the acting coil(s), the force experienced 

(i.e. , delivered) is not equal to that desired, i. e. , F.‘ j t F. . The specimen 

ii 

then does not experience the precise magnitude nor direction of the force 
required by Equation (1) and errors evolve. The magnitude of these errors 
and their effect on specimen motion have been simulated when we apply the 
command current to coil i, but the force delivered was computed on the basis 
of the off-axial forces due to the four-coil system, described in an appendix. 

2. 3. 2. 3 Analysis of Postulated System 

In terms of the block diagram (Figure 2. 3-1), when currents are calcu- 
lated at 4b or 4c, these currents must be used as input for electromagnetic 
field and force equations and the true resulting forces examined. Thus, the 
actual dynamics of the specimen via the spring equation and either the Fromm- 
Jehn equations (4b) or an approximation (4c) are studied relative to what is 
actually taking place. 

It will be shown that a workable, rather than an exact, solution is 
provided by the spring-on-axis approach. 


2-20 



The electromagnetic force acting on a sphere of radius on the axis 
of a coil carrying I amperes in N windings has been discussed previously* and 
is repeated here for convenience. 


F 



(NI) 2 A(y)G(x) 



y is the specimen-coil distance in units of the coil radius R^, x is the ratio of 
sphere radius to RF skin depth. The function G(x) has been discussed in 
Section 2. 1. 5 above. 


It has been found that an excellent approximation to the force near the 

2 

center of the tetrahedron is F = (NI) C 1 exp (C^) as shown graphically in 
Figure 2. 3-2 for the parameters: 


Coil Current 
Diameter of Sphere 
Material 

Distance between Coil Centers 

Distance between Center of Coil 
and Center of Tetrahedron 


45 amp (maximum acceptable) 
1 . 0 cm 
Al at 20°C 


3. 05 cm 
1. 87 cm 


Coil Diameter 2. 3 cm 

The constants C^ and C^ were determined by fitting the exact expression plotted 
in Figure 2. 3-2 over the control region near z = 1. 87. 


* Frost, R. T. , etal , Field Ma nagement for Positioning and Processing 
of Free Suspended Liquid Materials . General Electric Co. , Contract 
NAS 8-24683, Modification No. 2, Task IV, Final Report, May 15, 1970. 
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Axial Force, F - dynes 



z - cm 

Figure 2.3-2. Fromm-Jehn Force Equation and Approximation 




These conditions reproduce those for which more elaborate field calculations 
are available as well as laboratory measurements (see appendix for computer 
program). 

As indicated in Figure 2. 3-2, the approximation can be written: 

Log F = -1. 08z + 2, 93 (2) 

2 

dividing out (NI) : 

Log F/(NI) 2 = -l.08z-l.33 (3) 

or in terms of 4 (z = 4 + 1. 87): 

Log F/(NI) 2 = -1.084-3.36 (4) 

where 4 is referenced to the center of the tetrahedron. 

Consider the case of a single coil, given that a force F^ is desired along 
the axis. From Equation (4) with N = 3, I is easily calculated (and easily 
implemented electronically). Equation (4) agrees quite well over the control 
volume near z = 0 with the more exact solution for the single coil as shown in 
Figure 2. 3-3. 

Extending the approximation to include off-coil axis positions for two 
and three coils yields the comparison of Figures 2. 3-4 and 2. 3-5. The 
approximation to the left of zero (0) (toward the acting coils) fails completely. 
This is expected as it is in this region that the coil field superposition effects 
are most significant. On the positive side of zero (0), the approximation tends 
to agree rather well with the experimental data. 
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FORCE x 



DISTANCE (CM) 


FIGURE 2. 3-3 Force Along Axis of Symmetry of One Coil 


2-24 




DISTANCE 

FIGURE 2. 3-4 Force Along Axis of Symmetry for Two Coils Acting 
FORCE Simultaneously 



FIGURE 2. 3-5 


DISTANCE (CM) 

Force Along Axis of Symmetry 
Simultaneously 


tor Three Coils Acting 
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Note that although Equation (4) ove r-predicts forces to the left, forces 


are produced out to a distance of -1 cm, as in Figures 2* 3-4 and 2. 3-5* This 
is the basis of the 'workable 1 solution* Even though the precisely desired 
command forces are not delivered, the specimen does experience forces which 
tend to move it toward the origin* 

2* 3. 2* 4 Coil Vectors 

A symmetric orientation of the tetrahedron is shown in Figure 2*3-6* 

Let the edge length, A to B, be 2L. From the geometric properties of 
the regular tetrahedron, the unit vectors directed normally from the coils 
are : 

L 1 0 i - Jz/3 j* + y/if 3 k 
= -y/2/3 i + o j ->71/3 k 
L^ = y/2/ 3 i + o j - JTJT k 

L 4 = 0 f + ^jUl J* + Jl/3 k 

Edge AB I I to Y 
DC | | to X 

Face 1 = BDC = Coil 1 
Face 2 = ADB = Coil 2 
Face 3 = ABC = Coil 3 
Face 4 = ACD = Coil 4 

Figure 2. 3-6* Tetrahedron Orientation 
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Assume each coil is capable of producing force F. on a particle, at 
rest, at vector position p{x,y > z). It is desired to determine the force F^, 
which each coil must produce so as to achieve a net force in the desired direction, 
i.e. , in the direction opposite to the position error p. The physical problem is 
to return the particle to the origin (0, 0, 0) by manipulation of the coil forces. 

The direction taken by the particle is determined by the ratio of the coil forces. 

If we define a unit vector oppositely directly to the position error, 

p # - -2- the vector problem is then to solve for F. from the set of equations 
P 1 

p' = L F. L. , i = 1, . . , 4 (5) 

.ii 

i 

The direction cosines of the unit vector p are cos O', cos cos y . We require 
then that the coefficients of i , j , k in p 7 be equal to the sum of the i, j , k, 
coefficients of the L' s. 

Let a = y/Z/3 , b = V 1 / 3 

Then 


COS O' 

i = 

^(F 1 (0) ) + (-a F 2 ) ■+ (a F 3 ) + (0 F 4 ) ^ i 

(6) 

cos /3 

j = 

^■aFj) + (0F 2 ) +(0F 3 ) +(aF 4 )Jj 

(7) 

cos y 

k = 

f(b F l ) + (-b F z ) +(-bF 3 ) +(bF 4 )'jk 

(8) 


The vectors L. are linearly dependent, but any combination of three of 
them are linearly independent. Consider an arbitrary value for p . In the 
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physical geometry, p' can be represented uniquely by selecting the three 
appropriate coils of the four. That is, a maximum of three of the pushing 
(positive F^) coils will return the particle to the origin. The fourth coil is not 
employed. Equations (6), (7) and (8) will form a 3 x 3 system of equations if 
one of the F. 1 s is set to zero. 

l 

The system in Equations (6), (7) and (8) are written in standard matrix 
form as 



(A 1 ) 


Since these three equations do not yet determine the four force components F. 
one more relation involving the F' s is required. We may excite the coils in such 
a manner as to give a fixed arithmetic sum for the four force components: 



( 10 ) 


( 11 ) 
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Set C to 1 . Then 


a o 


cos OL 


o a 


cos /S 


-b a 


1 1 


cos y 


F = A -1 p' 


Substituting a and b in A and inverting yields 


Where 


c = 0.61237273 


d = 0.43301290 


e = 0.25000000 


Thus, given p'(x, y, z), F is determined from Equation (13). 
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A number of solutions to equation (12) have been obtained and can be 


summarized as follows: 


Case 

For p' 

General Result 

(1) 

Not in the plane of 

one - F force, three + F forces 


two coils 


( 2 ) 

In the plane of two coils 

two equal - F forces, two + F 



forces 

( 3 ) 

Along one coil axis 

three zero forces, one + F force 


The negative F forces occurring in Cases (1) and (2) are due to the fact 
that system of Equation (13) is dependent and there is no way to mathematically 
restrict the signs of F.„ The negative forces are very easily handled. Since 
any one of the F 1 s can be resolved into components on the other three, the 
negative F of Case (1) can be resolved onto the three positive F. 1 s. It can be 
shown that this resolution is accomplished by simply adding directly the 
absolute magnitude of the negative F to the magnitudes of the other forces. 

Note than when this manipulation of the negatives is performed, the constraint 
Equation (11) is broken; i.e. , 2 F. ^ C or 1 as for Equation (12). This is un- 
important as it is the F ratios that are desired at this stage. 


F P 

The above analysis is the "coil solution" to the equation — = -k — , 

where k = 1. The second term (damping term) of the equation (f 


dt 


velocity = 

is handled in an analogous manner. The velocity data provided by the 


position sensor electronics, x, y, z is resolved onto the coil-origin axis and 
the vector sum of the position term and velocity term is taken along each axes, 
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*' e '’ ^i(position) "^i (velocity) 

Again, any negative forces arising are resolved onto the positive (pushing) 
forces. 

2. 3. 2. 5 Resulting Control Equations 

The control equations define the current, I , to be delivered to each coil. 
To perform this calculation we must first calculate the negative of the position 
vector p resolved on each coil vector, A.. These are given below for the 
tetrahedral geometry. 

A = + ^2/3 y - -s/l / 3 z 

A = V 2/3 x + n/ 1 / 3 z 

A^ = - \JZl 3 x + Vl /3 z 

A^ = - V 2/3 y - n/ 1 / 3 z 

These functions are differentiated to obtain A^ , for the damping term. 

The force required of each coil is now calculated from the above using the 
spring-mass -dashpot form of solution 

* 

F. = +K A. + b A. 

l ii 

As discussed previously, the magnitude of negative forces are added to each 
of the other three forces so that a positive (pushing) solution is obtained. The 
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force required is at this point to be considered a constant and the current for 
each coil is given by the approximation to the Fromm-Jehn solution by 

(m + n) 

I. = F. e v 1 
1 1 

where m and n are scaling constants. These constants are chosen so that 
maximum current is delivered before a specimen approaches the minimum 
force position among three coils. 

2. 3. 2. 6 Mechanizing Electronics 

A complete electronic block diagram and schematic is contained in 
Appendix A. The implementing electronics of the coordinate conversion 
equations is depicted in Zones B10-D10 to B8-D8 in GE Drawing ER 47J225305, 
Sheet 1 in that appendix. The Force /Current Computation electronics are 
shown in Zones A8-D8 through A3-D3 on the same drawing. As noted in that 
appendix, the design effort on these circuits had not been completed nor test 
circuits built at the time of contract termination. 

2. 3. 3 Analysis of Motion of Aluminum Specimen in Four-Coil System 

A limited number of computer simulations were made to study the capture 
capability of the four-coil system. Each coil was limited to a maximum 45 
amps. The analytic field solution was employed in the analysis. Although 
differences exist between these solutions and the experimental measurements, 
the two agree well enough (see Figures 2. 3-3, 2. 3-4, 2. 3-5 to warrant use of the 
field solution. 
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The specimen material chosen was an aluminum sphere of 1 cm diameter. 
An arbitrary temperature of 700 C was assigned and the skin depth calculated 

at this temperature is 0. 04 < 6 < 0. 05 cm depending on coil frequency. The 

5 5 5 5 

frequencies assigned were 10 , 1.16 x 10 , 0.78 x 10 , and 0.84 x 10 Hz. 

2. 3. 3. 1 Sting Position 1 

The specimen position selected for the first runs was that shown in 
Figure 2.3-7 with the coil number designations as indicated in Figure 2.3-6 
This was one postulated position in effect when work on the contract was 
terminated. However, this position is unacceptable since the action of the 
postulated control is to eject the specimen completely from the control volume 



Figure 2. 3-7. Sting Position 1 
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At the position location of the specimen as given above, off-axis effects 
of Coil 3 are significant. Presuming that the direction of first motion of the 
specimen is in the direction of the total force vector given by Gray/Wouch, 
the specimen was ejected from the system for initial velocities of 0„ 0 and 0.5 
cm/ sec. 

This result was due primarily to the force component perpendicular to 
the Coil 3 center line. The initial motion is indicated by V in Figure 2.3-7. 
Note that Coil 1 cannot aid in establishing motion toward the origin (0, 0, 0). 

The specimen is thus moving in the volume defined by Coils 3, 2 and 4 (Coil 2 
and Coil 4 are not shown). This combination puts the specimen in the weakest 
position of the system, i. e., approaching the three coil symmetry line where 
"negative" (expulsive) forces exist beyond about 0.65 cm (Gray/Wouch results) 
or 1 cm (lab results). The specimen eventually enters the "negative" force 
volume and is expelled. 

2. 3. 3. 2 Sting Position 2 

A second specimen position was chosen so that the axis of the sting 
lay along the Coil 3 center line and was 1 cm from the origin (0. 87 from Coil 3). 
The initial motion here is along the center line. 

Results of several runs from this position are shown in Figure 2. 3-8. 

Up to 3 cm/ sec. initial velocities are managed by the system such that critical 
damping is achieved. Between 3 and 3.5 cm/sec, the specimen overshoots the 
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Specimen 0.87 cm in front and on axis of Coil 3 
Aluminum Specimen; 1 cm dia @700 C 
Coil-to Coil distance = 3.05 cm 
V Q directed toward center (0, 0, 0) 

Maximum current per coil = 45 amps 



Figure 2.3-8. Motion from Sting Position 











origin but is eventually turned around and damped. Between 3.5 and 4. 0 
cm/ sec, a velocity exists such that the three damping coils (Coils 1, 2 and 
4) cannot handle the initial velocity if limited to 45 amps/coil. The result 
is ejection of the specimen, 

2. 3 0 3 0 3 Sting Position 3 

The unacceptability of Sting Position 1 indicated that a change in the 
sting location and orientation was necessary. Sting position 2 yielded well 
controlled specimen motions for reasonable initial velocities. No hard data 
were available on "wetting" properties of various candidate materials and 
therefore initial expected velocities could not be accurately predicted. The 
sting position chosen for this simulation is sketched in Figure 2. 3-9 and 
represents a mechanically feasible configuration as compared to Sting 
position 2. The sting enters the control volume between Coils 2 and 4 (not 
shown in Figure 2. 3-10). This location would permit a relatively strong 
force to be developed by Coils 2 and 4 working in concert to lift the material 
from the sting in case it "sticks" due to some partial wetting of the sting 
platform by the molten specimen. The initial velocity is, in effect, aimed 
toward the center of the control volume to permit effective utilization of 
Coils 1 and 3 as retarders. 

2. 3. 3.4 Ejection Geometry and Constraints for Sting Position 3 

1. Coil to coil separation: 3.05 cm 

2. Sphere radius: 0.5 cm 

Coil radius: 1. 15 cm 
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Figure 2.3-10. Force Function for Sting Position 2 



3. Release: sphere off parallel to longitudinal axis of sting and toward 

origin. 

4. Coil forces: Data are shown in Figure 2. 3-10 for two coils with 
I = 45 amps, N = 3. 

5. Problem geometry: Figure 2. 3-9. 

6. Material properties (molten): 


t 

>, Skin Depth, cm 
at 10 5 Hz. 

Mass, gm 

(1cm dia. spheres) G (x) * 

Silver (Ag) 

0. 066 

4. 97 

0. 802 

Copper (Cu) 

0.073 

4. 29 

0. 781 

Aluminum (Al) 

0. 078 

1. 24 

0. 901 

Iron (Fe) 

0. 188 

3. 58 

0. 428 

Nickel (Ni) 

0. 147 

4. 48 

0. 555 


*See Section 2. 1. 4 
2. 3. 3 . 5 Procedure 

From Figure 2. 3-9, the kinetic energy imparted to the sphere at the sting, 
2 

namely, 1/2 mv , must be removed by the force experienced by the sphere over 
the distance 2(1. 08). This is equal to the average force over the distance times 

the total distance, Al. 


0.5 mv 2 = J F(l) dl 

/ f 

(16) 

- -J l 

/ -l J 


= G(x) I 2 f 
L -i 

C(l) dl ^ \ i 1 dl Al 

(17) 
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The bracketed term on the right can be evaluated from Figure 2. 3-10 by 

numerical integration of F (from -1, 08 to 1. 08) and then dividing out G(x) and 
2 

I „ The average force exerted is 

F t (1) = G(x) I 2 3, 48 x 10' 3 (18) 

Here mixed units are used, the force in dynes and the current in amperes. 

In Figure 2„3-10, F ~ 6.4 dynes. Equation (18) is used instead of the 
more familiar*: 

F = 3/50 (NI) 2 A(y) G(x)(R2/Rl) 3 in MKS units. ( 19 ) 

Figure 2. 3-10 is presumed to be the more accurate description of the force 
function due to the two "stopping" coils, 1 and 3. 


From Equations (16) and (17) 


V 

o 


2. *3. 48 x 10" 3 * 2. 16 * 


G(x) 


m 


1/2 


and 


1/2 

I = 8. 16 (m/G(x)) W Y 

o 


( 20 ) 


( 21 ) 


* Frost, R. T., et al, Investigation of the Preparation of Materials in Space, 
Task IV, Field Managem ent _for Positioning and Processing of Free Suspended 
Liquid Materials . Final Report (U), GE-SSL for NASA, Contract No, 

DCN1 -9-54-20055, 52, 15 May 1970, Pg. VII-5. 
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2. 3. 3. 6 Results 


Figure 2. 3-11 shows the application of Equation (21) for the five molten 
materials examined. 

As may have been expected, aluminum, because of its light weight and 
small skin depth, is the least current demanding. At any ejection velocity V , 
the other four materials (Cu, Ag, Fe and Ni) require from 2 to 2-1/2 times 
the current necessary to stop an A1 sphere. 


For a maximum coil current of 45 amps, the maximum tolerable ejection 
velocity for the five materials are 


Maximum Sting Velocity for 

Molten Material Maximum Coil Current 45 amps 

(1 cm sphere dia) (cm/ sec. ) 


A1 


4. 37 


Cu 


2. 36 


Ag 2. 22 

Fe and Ni 1 . 92 


Note that the ejection geometry as postulated confines the sphere to the 
plane containing the axes of the two stopping coils. Should the sphere have a 
velocity component out of this plane, Coils 2 or 4 (not seen in Figure 2. 3-9) 
will act (via the servo system) to restore the sphere to the plane. This force 
will however tend to increase the apparent sphere ejection velocity, "fighting" 
Coils 1 and 3. If ejection parallel (or nearly so) to the sting longitudinal axis 
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cannot be assured, the results herein are to be viewed on the optimistic side. 
That is, the in-plane case represents the best (current conservative) configura- 
tion for handling an ejection velocity, v. 

2. 3. 4 Specimen Capture 

After each specimen has been melted, and then has again become solid, 
it must be removed from the positioning coil system to make way for the next 
specimen to be melted. To do this, the position control system places the 
specimen to be ejected in a location from which one or more coils acting 
simultaneously can push the specimen from the interior of the set of coils. 
Several techniques were considered to prevent ejected specimens from returning 
to the vicinity of the coils and electron beam once they had been ejected. They 
are listed below using the words capture or non-capture depending upon the 
degree of confinement resulting from each technique. 

2.3.4. 1 Non-capture, Freely Drifting 

During the course of the experiment the entire apparatus will experience 
an acceleration due to the drag of the atmosphere upon the vehicle, due to 
gravity gradient effects and due to rotation of the vehicle (if any rotation exists). 
Any acceleration relatively constant in direction will cause freely drifting 
specimens to drift to one side of the vacuum chamber and to congregate there. 

If the net acceleration upon the experiment is in a suitable direction for a 
sufficiently long time the ejected specimens would remain away from the coils 
and electron beam. 
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2„3. 4.2 Non-capture, Inelastic Rebound 

The technique described above could be improved by decreasing the length 
of time it takes each specimen to attain a relatively motionless state if an energy 
absorbing surface were placed in the path of the ejected specimens to deflect and 
reduce their velocity or if an energy absorbing surface were placed somewhere 
in the vicinity of the region in which they congregate. This latter surface might 
be a sticky surface of a suitable material, such as a low vapor pressure vacuum 
grease, to be certain of each remaining in that region. Materials used to 
decrease the speed of ejected specimens would have to tolerate contact with the 
hot metal spheres and the degree of outgassing produced by this contact would 
have to be evaluated to determine to what extent such outgassing would necessitate 
delaying the melting of subsequent specimens. 

2. 3. 4.3 Capture, Fiber Trap 

The use of very fine metallic fibers to construct a valve-like device was 
considered at length in response to the desire for a simple passive means of 
capturing processed specimens on or in the wheel to which they were initially 
attached. In this technique one or more fibers are installed around the mouth 
of a cavity. If a specimen ejected from the coil system possesses sufficient 
kinetic energy to deflect and pass the fibers but insufficient energy to rebound 
out, each specimen could be stored in the wheel in the same order in which they 
were initially while yet attached to the wheel. 
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Several tests were performed to assess characteristics of such a tech- 
nique by using a pendulum-like apparatus which permitted a metal sphere of a 
size comparable to that used for the specimen wheel study to swing into a 
"capture port" at several known speeds,, The result of these tests indicated 
that kinetic energy at least of the order of 40 ergs is required to operate such 
a capture device successfully. This concept is therefore marginal for use 
with the low forces dictated by the very limited power availability for this 
experiment from the M512 facility. 

In order to carry a greater number of specimens on each wheel the melted 
specimens might be collected in a separate container thereby allowing the space 
formerly used for stowage cavities to be used to attach more specimens to the 
wheel. The fiber trap technique could also be applied tp such a container. 
However, because it does not assure passage through the port in the event that 
one or more previously captured specimens are in or near the port, thereby 
blocking the entry of another specimen, it was abandoned for this application 
in favor of one of the following techniques. 

2. 3. 4. 4 Capture, Coil 

A desirable technique for operating a port in a collection box is to place 
a coil about the port and operating it, in a manner similar to that used for the 
positioning coils, to repel previously captured specimens from the port. The 
coil could be turned on only as commanded by "specimen- sensor s" (e.g. , light 
emitting diodes paired with light detecting diodes) or it could be on continuously 
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to be turned off by logic circuitry whenever a new specimen is to be sent into 
the box. Such a coil would have to be placed far enough from the positioning 
coils to prevent the mutual induction between it and any one of the positioning 
coils from being large enough to cause a disturbance in the servo-positioning 
system. The passive techniques, described in Sections 2. 3. 4. 1 and 2.3.4. 2, 
could be applied to the captured specimens to minimize the chance that an 
incoming specimen could collide with a previously captured specimen in such 
a manner that it could rebound out through the port before the coil could force 
it into the container. 

2. 3. 4. 5 Capture, Door 

An alternative to the preceding capture, coil technique involves the sub- 
stitution of a door or gate which opens inward to accept new specimens at the 
same time pushing previously captured specimens away from the port. 
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2.4 FOUR-COIL. ENGINEERING DEVELOPMENT BREADBOARD 


In order to perform functional tests of the four-coil configuration involving 
forces, coil interactions, position sensing (electromagnetic, infrared and 
optical), impedance matching circuits and general servo behavior, a bread- 
board was constructed. The construction was done simultaneously with other 
studies such as position sensing, coil optimization, mechanical considerations, 
etc. The most important contribution made by this breadboard was a full three - 
axis control demonstration of a free and neutrally buoyant hollow aluminum 
sphere. Demonstrations were performed and a movie was taken (included as 
a part of this effort) of this sphere immersed in a mixture of silicone oils or 
in a brine solution. The set-up can be seen in the photograph, Figure 2.4-1, 
where the white ball (representing the levitated specimen) is visible centered 
in the liquid filled glass sphere. 



Figure 2.4-1. Engineering Breadboard Coil and Specimen Arrangement 
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2„4„1 Electronic Circuitry 


The electronics conformed closely to the circuitry for the proposed 
system in all respects except for the infrared position sensing and the 
coordinate converter. These last items were not complete at that time so they 
were not incorporated. Instead, the ball corresponding to the specimen was 
painted white and illuminated by several incandescent lamps. Position sensing 
was accomplished by three single axis sensors which required no coordinate 
converter. The sensors and incandescent lamps can be seen by careful examin- 
ation of Figure 2.4-1. 


2.4.2 Position Sensing 

The single axis sensor looked into the enclosure between two coils at their 
closest approach (mid-way between vertices) and was aligned so that the 
position error was sensed along a plane containing the sensor and the centers 
of the two opposite coils, as shown in Figure 2.4-2. 


Sensor 




Figure 2.4-2. Direction of Sensed Error on Breadboard 
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In Figure 2.4-2 the arrow shows the direction of the sensed error. This error 
drives either Coil 1 or Coil 2 depending on whether the specimen is above or 
below the center of the facility. 

There are two other sensors similarly placed between other pairs of 
coils so that position information is available on the specimen throughout the 
coil volume for a two centimeter radial displacement of the specimen from the 
center of the coil facility. 

If the preamp design for infrared position sensing had been completed, 
the four -coil breadboard could have been reworked to use it with the single 
axis sensors, both the single axis and dual axis (4 quadrant) sensors being 
silicon junction devices. To use infrared position sensors, a stainless steel 
bob was to have been suspended on a long pendulum and heated in an induction 
heater. The hot pendulum bob would then be transferred to the four-coil 
facility for tests. 

2.4. 3 Test Procedures 

All test runs using the four-coil breadboard were with the neutrally 
buoyant sphere. Four switches were added to the position sensing circuits to 
displace the ball in the facility which can be seen at the upper right hand 
corner of Figure 2.4-3. Actuating a switch caused the sphere to move in an 
arbitrary direction. Releasing the switch enabled the servo to re-center the 
sphere in the coil facility. Four switches provided four directions for sphere 
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Reproduced from 
best available 


copy. 



Figure 2.4-3. Complete Four-Coil Engineering Breadboard 
displacement. This displacement was accomplished by introducing a position 
voltage in one of the position sensor amplifiers for the servo to null out by 
moving the sphere away from center. Simultaneously, the other position 
sensor amplifiers are disabled. 

2.4.4 Schematic 

A full schematic is given in GE Drawing ER 47E225310 included as 
Figure 2.4-4 in this report. Figure 2.4-5 depicts the wiring for the positioning 
switches and is used with Figure 2.4-4. 
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2. 5 POSITION SENSING AND SERVO SYSTEM 


2, 5. 1 Optical Position Sensing 

The location of the hot specimen inside the coil system was determined, 
in the reference design, by a four quadrant sensor in a suitable optical arrange- 
ment that detected energy radiated by the specimen in the near-infrared spec- 
trum,, Electronic preamplifiers and a divider, to normalize the amplitude 
variation of the cooling ball, provided signals to the coordinate converter 
proportional to errors along three orthogonal planes. 

2. 5„ 1„ 1 Molecular Environment 

The position sensing system was to operate in a vacuum chamber. Since 
the experiment consisted of melting metal spheres, metal vapor would be the 
most significant contaminent during the experimental period. Calculations 
were made to determine the thickness and optical properties of the metal 
films deposited on surfaces in such an environment. For aluminum spheres, 
a metal film of between 10 angstroms and 100 angstroms thickness would be 
deposited on exposed hardware surfaces during a period of approximately one 
second at a radial distance of two centimeters from the source. The melting 
phase of the experiment was expected to last for perhaps two seconds for each 
specimen, and since metal films of 10 angstroms to 100 angstroms thickness 
are opaque to the visible and near infrared radiation, it was evident that trans- 
missive optical components, if used, could not be directly exposed to the metal 
vapors. 
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2. 5. 1. 2 Thermal Environment 


The electron beam used to melt specimens would be operated at a power 
of about 1. 6 kilowatts for between two and four seconds per specimen. During 
this period of time, between 4000 and 8000 joules of energy would be deposited 
in the specimens and the attached sting. This energy must be either radiated 
or conducted away by equilibrating physical processes. The nearness of the 
optical system and its associated housing to the target specimen would influence 
the portion of the total energy that would be absorbed by these components. A 
large portion of the energy absorbed in the specimen would be conducted away 
through the sting to the relatively massive support wheel and associated support 
mechanisms (see Figure 1.1-1). The remaining mode of energy transfer would 
be through a radiative process. Fortunately, the thin layer of metal vapor 
deposited on the optical system and associated housing will effectively reflect 
the incident thermal radiation falling on the sensor subsystem. Since, how- 
ever, all the other parts of the chamber undergo similar vapor treatments, the 
amount of radiation which would be received by the optical subsystem was 
difficult to evaluate exactly. 

Calculations were made in which it was assumed that all the energy 
radiated from the hot spheres into a solid angle subtended by the proposed 
optical systems were absorbed. The optical subsystems were in good thermal 
contact with a large heat sink. The calculation showed that we could expect 
only a few degrees centigrade rise in temperature within the optical housing. 
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Therefore as long as reasonable precautions were taken to ensure a good heat 
sink for the housing of the optical components and that the housing itself be 

V 

made of highly conductive material there would be little problem with thermal 
conditions in the optical positioning devices, 

2. 5, 1, 3 Thermal Gradients 

Both the positioning sensors and the temperature measuring sensors 
(Section 2, 6) would be affected by thermal gradients within the specimen. In 
a NASA film showing the melting of a 0, 5 cm diameter sphere, temperature 
gradients within the sphere were clearly visible. This could result in a false 
interpretation of the position of a specimen. A one centimeter diameter 
specimen would perhaps increase the magnitude and spatial extent of the thermal 
gradients owing to the eight-fold increase in its mass and a decrease in the 
ratio of its area to volume. The control system as presently conceived would 
force the specimen into a position which corresponds to the centroid of the 
energy emitted from the specimen. 

2. 5. 1.4 Optical Line -of-Sight 

A constraint on the optical line -of-sight for the position sensing system 
was that within the tetrahedral geometry of the coils there be a spherical 
control volume of four centimeters diameter. If the specimen was anywhere 
within this control volume the position sensor must not have an ambiguity in 
its control function output. Stated simply this means that the specimen must 
be controlled back to the midpoint of the volume and not ejected if it was 
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contained within this volume. A second condition which must be taken into 


account was that the housing for the optical system should not be closer than 
approximately two coil diameters from any particular coil so as not to in- 
fluence the field about this coil. In addition we desired the fewest required 
number of sensors so as to minimize the number of components of the system. 

Figure 2.5-1 shows the schematic configuration of the location of the two 
optical positioning sensors. Each sensor consisted of a pinhole lens system, a 
folding mirror, a detector, and optical filter. These sensors were placed at 
two of the corners of a tetrahedron configuration; the control coils of which 
were on the faces of the tetrahedron. The size of the detector and the "focal 
length" of the pinhole lens was chosen so that an inscribed sphere of four 
centimeters is completely within the optical field-of- view of the sensor and 
there are no areas which vignette within this field-of-view„ The two sensors 
were used to define a three-axis cartesian coordinate system. The detector 
within each sensor was a four quadrant PIN-Si photodiode. A four quadrant 
detector can be used to define a two-axis coordinate system within the plane 
of the detector by judiciously summing the outputs in the proper fashion from 
each quadrant of the detector. When the image of a sphere is exactly on the 
center the outputs from all four quadrants are equal and a null signal can be 
obtained. By employing two such sensors we can define the three -axis 
coordinate system within the tetrahedral geometry as follows. 
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Figure 2.5-1. Optical Sensor Orientation 


Figure 2. 5-2 is the same as Figure 2. 3-6 but is repeated here to 
simplify the explanation of the relation of the optical sensors to the coordinate 
system utilized. 


Edge Length = 2L = S 
Edge AB || Y 
DC || X 

Face 1 = BDC = Coil 1 
Face 2 = ADB = Coil 2 
Face 3 = ABC = Coil 3 
Face 4 = ACD = Coil 4 


Figure 2. 5-2. Tetrahedron Orientation 
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Locating an optical sensor at point B (O, L, L ) we can rotate the detector 

so that the X axis, that is a local coordinate axis within the detector, is 
B 

parallel to the X axis within the tetrahedron. Any motion of the specimen 
along the X axis within the system results in movement of the specimen image 
along the axis of the detector in accordance with the following equation: 

X = (M/f)X B . 

Where M is the distance from the pinhole lens to the center of the coordinate 
system; f is the "focal length" of the pinhole lens and equals the distance from 
that pinhole to the detector; and X Q is the position along the axis in the 
detector plane parallel to the tetrahedron X axis. Another detector located at 
point C can be oriented so that its axis is parallel to the Y axis within the 
tetrahedral geometry. A similar equation is derived for the motion of the 
image on this detector in terms of motion within the tetrahedron; i. e., 

Y = (M/f)Y c . 


Motion of the specimen along the Z axis within the tetrahedron is given by a 
combination of the motion of the image on the detector at point B and the per- 
pendicular axis (which axis lies in the plane defined by the system X and Z 
axes) within the detector at point C, and the following equation is derived. 


Z = 


Z c (M/f)/Sin 54°44 


Z 

c 

f 


Cos 54°44 


,4. , 

Cos 54°44 + Sin 54°44 


K] 

r> 1 n ' 


.Sin 54 44 - —h Cos 54 44 J 



• f 
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The electrical signal for each axis of motion is determined by the responsivity 
of the detector to motion of the image on the sensitive element. 

Figure 2. 5-3 shows the relative output from the detector packages versus 
displacement of a specimen from the center of the tetrahedron. This function 
is single valued for all points within the control volume for a given specimen 
temperature. Since the temperature of the specimens would vary, we must 
normalize the output as a function of the absolute output from the detectors as 
will be discussed in Section 2.5. 1.7 below. 

2. 5. 1. 5 Signal-to-Noise Ratio for Position Sensors 

Since multi-element arrays were necessary to define the coordinate 
system within the field-of- vie w and, in the early considerations of detectors 
we were also attempting to measure temperatures of the specimens with the 
same detectors, several different detector types were considered. A multi- 
element array of detectors of perhaps PbS and PbSe were found to be both 
expensive and long lead-time items. The temperature sensing element was 
therefore separated from the position sensing element employing the relatively 
low cost off-the-shelf silicon multi-element arrays for position sensing. 
Temperature sensing development is discussed in Section 2. 6 below. 

The signal from a silicon detector placed at the foregoing positions in the 
tetrahedral geometry is controlled by the absolute temperature of the specimen, 
the specimen emissivity, the aperture of the pinhole lens, and the responsivity 
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of the detector element. For the purpose of calculation, the average 
emissivity from typical candidate specimens was assumed to be approximately 
0. 1. The size of each of the specimens was to be one centimeter in diameter. 
The aperture of the lens was assumed to be two millimeters and the responsivity 
of each of the detectors from published data by United Detector Technology is 
0.5 amps per watt. All of the candidate specimens were estimated to solidify 
by about 800°K and this defined the lower temperature limit for the control 
system. Conversely all candidate specimens would melt below a temperature 
of 177 3°K defining the upper temperature limit and demarcating the dynamic 
range requirements of the position and temperature sensors. The calculations 
were performed for a specimen at the null position in which the sum of two 
detectors each seeing one quarter of the specimen were added tbgether. A 
minimum output current was calculated to be 3. 3 x 10 ^ amps. 

The noise associated with each of the detectors is controlled by the 

electrical bandpass of the system and the leakage current at the silicon junction. 

A safety factor of ten from the published data was incorporated into the design 

figures. The output electrical noise in a 100 Hz bandpass is on the order of 

8 x 10 ^ amps. An additional contribution to noise would be false signals 

from elements within the field-of-view not associated with the specimen. 

For calculation of these signals it was assumed that the background consisted of 

a blackbody at 300 K with an emissivity of unity. The output resulting from this 

- 1 9 

assumed background would be 4x10 amp and therefore could be 
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totally neglected. Figure 2.5-4 shows the signal-to-noise ratio for two cells 
versus the temperature of candidate specimens for the null position. It is 
worth noting that at a temperature of 800 K, a signal-to-noise of approximately 
ten was expected. Since in this application the detectors were operating in a 
differential mode, the accuracy of the null position would be controlled by the 
noise output from the detector elements. Attention must be called to the very 
low current levels expected and to the extreme precautions required in the 
design of the electronic shields. 

A definite improvement in signal-to-noise ratio could be achieved through 
the use of a mechanical chopper, such as the Bullova Tuning Fork Choppers. 
However, questions of reliability in the functioning of a chopper in this 
metallic environment led us to attempt to operate in a dc fashion. The design 
however does not preclude the use of such a chopping system and if actual 
testing of the breadboard equipment would demonstrate the necessity of im- 
proving signal-to-noise ratio then such a chopper could be added at a later 
time. 

2. 5. 1. 6 Mechanical Design 

An engineering model of a combined position and temperature sensor 
was constructed and housed in an aluminum case which is 1. 15 by 4. 5 inches 
long (exclusive of connectors) by 1. 2 inches high. Figure 2. 5-5 is an outline 
drawing of the prototype assembly and Figure 2.5-6 is a photograph of the 
unit fabricated together with the temperature sensor electronics in breadboard 
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Temperature - °K 

Figure 2.5-4. Signal-to-Noise Ratio of Silicon Detector Elements 
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Figure 2.5-5. Position and Temperature Sensor Housing 



Figure 2.5-6. Prototype Sensor and Electronics Breadboard 
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form. Each sensor, that is the position and temperature sensor, sees the 
entire control volume by means of a pyramidal mirror structure and though 
its own aperture. The mirrors are utilized to fold the optical path and permit 
both sensors to view the control volume. They also serve to introduce a re- 
flective optical element in the optical path from the hot, potentially evaporating 
specimen. Any vaporization deposits will have a negligible effect upon the 
performance of the front-surface mirrors and because of the folded optical 
paths, they would not reduce the transmittance of the refractive optics. 

One aperture serves as the pinhole lens for the position sensor, the second 
aperture simply limits the total field of the temperature sensor to be some- 
what larger than the control volume. An optical filter in each path rejects 
the visual portion of the optical spectrum. This filter is used to prevent 
spurious signals from entering the detectors due to the photo floodlamp used 
in conjunction with the motion pictures that were to be taken during the 
melting experiments. At each end of the case is a connector which connects 
the respective detectors to the appropriate electronic circuitry. The 
mechanical design serves as almost a continuous electro-magnetic shield to 
prevent spurious radiated signals from entering the signal channels. 

Only one of the sensor elements would contain both temperature and 
position sensing elements. The second sensor would contain the position 
sensing element alone. As noted in Section 2.5. 1.4 above, the other sensor 
is required to supply the third axis of positioning information. 
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2. 5. 1. 7 Electronic Circuitry Design 

An optical schematic of the quadrant detector and aperture, without the 
flat mirror, is included in Figure 2. 5-7. 


+Y 



Figure 2.5-7. Optical Schematic for Y and Z Axes 

The specimen is shown displaced in "Y" so that detector quadrant -Y is 
illuminated more than +Y. The photo- generated current in -Y exceeds that 
in +Y by the amount of the error. Note that all the radiation passing through 
the aperture falls on the four quadrants. This requirement must be met in 
order to normalize the amplitude variation due to the cooling specimen. 
Normalizing for changes in temperature and emissivity is done by dividing 
the differential output of opposite quadrants by the sum of all four quadrant 
outputs. This scheme is depicted schematically in Figure 2.5-8 and the 
proposed implementation is shown in Zones B10-D10 to B12-D12 on Drawing 
ER 47J225305 of Appendix A. As noted in the preface to that appendix, these 
circuits had not been completely developed. 
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Figure 2. 5-8. Conceptual Schematic of Normalizing Circuits 

The output of the dividers is accurately proportional to the specimen 
displacement so long as the specimen temperature is high enough to produce 
a workable signal, and all radiation through the aperture falls on the detector. 

The arrangement for the X axis is different in that two quadrants are 
parallel as shown in Figure 2. 5-9 so that +X is from these quadrants and -X 
is from the opposite two quadrants. Implementation of these circuits is 
included in the Appendix A in the same area as mentioned immediately above. 

Laboratory measurements using a blackbody source (IR Industries 
Model 406) produced a preamp output of 10 volts at a source temperature of 
1000°C and 0. 035 volts at 590°C. 
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Figure 2. 5-9. Optical Schematic for X Axis 

2. 5. 2 Electromagnetic Position Sensing 
2„ 5. 2. 1 Basic Concept 

The presence of the specimen in the alternating magnetic field gives 
rise to two changes in the coil' s properties. The specimen will absorb power 
from the field and, due to the induced magnetic dipole, it will decrease the 
inductance of the coil. Examination of Figure 2.5-10 shows that both these 
factors tend to reduce the power factor angle 0 of the coil. Both the above 
effects are increased as the field amplitude increases. Consequently, if 
the specimen approaches a coil the power factor of that coil will be reduced. 
The power factor angle then could be used, under the proper conditions, 
to determine relative ball position within the coil system. 

The electromagnetic position sensing system developed under this 
contract sensed the power factor angle of one coil and compared it to the 
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Figure 2. 5-10. 



average of the other three coils to define the specimen position relative to the 
coil system. The block diagram of this system is shown in Figure 2.5-11. 

Figure 2. 5-12 depicts the waveforms and the reference pulse which 
represents 0. The dotted curves show the waveforms when the specimen is 
close to the coil. 

2. 5. 2. 3 Operation 

The voltage and current waveforms are normalized by using a zero 
cross-over circuit providing from that point on the bipolar switching signals, 
preserving only the phase of the two waveforms. V and I are defined as the 
positive half-cycles of the waveforms and V and I are the negative half-cycles 
of the waveforms. The logic provides a reference pulse only when the two 
waveforms are of opposite polarity. All reference pulses from all coils are 
of the same height and are exactly proportional to the respective phase angles 
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in width. If all pulses have the same average value, then all coils have the 
same power factor angle and the specimen is centered. Under this condition, 
the position error to the servo of each coil is zero because 


'K0, 


K0. 


K0 


K V 


+ 


= 0 = error of coil #1. 


If the specimen moves toward coil #1, then K0^ is reduced and K©^, 
K0^ and K0^ are increased and 

K0 K0 K0 ' 

— + — + — ^ < °* 

so coil #1 has a negative position error. At the same time coils 2, 3 and 4 
have a positive position error. In the servo, a negative error increases the 
power amplifier drive so as to repel the specimen from that coil. 



2. 5. 2. 2 Electronic Circuitry 

The basic scheme of electromagnetic position sensing had been devised 
and demonstrated under Contract NAS 8-24683. Under the aegis of the present 
contract, the concept had been developed for use with the tetrahedral coil 
configuration. The schematic of a single channel the circuit developed is 
depicted on Figure 2.5-13. A breadboard incorporating electromagnetic 
sensing channels for all four of the tetrahedral coils was fabricated and is 
shown in Figure 2.5-14. Check-out was completed for the channels associated 
with the three coils in the breadboard which had upward directed force 
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Figure 2. 5-14 Breadboard of Electromagnetic Position Sensing Circuits 



components before contract termination. With these three channels operating, 
the capability of the system was demonstrated to move a simulated specimen 
suspended from a fine wire. No effort was exerted on the fourth channel 
since it was not utilizing in controlling a specimen suspended as a pendulum 
in the laboratory. 
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2.6 TWO-COLOR PYROMETER 

2. 6. 1 Overall System Description 

Calculations were performed to show the feasibility of using a dual 
detector system consisting of a silicon detector with a lead sulphide cell 
immediately behind it. This detector system was purchased from Electro 
Nuclear Laboratories (ENL), Menlo Park, California, and made into a sensor 
head package that included a small circular aperture, deflecting mirror and 
IR filter as sketched in Figure 2.6-1. 


Hot Specimen 


Aperture^" 



Figure 2.6-1. Two-Color Pyrometer Block Diagram 


The radiant power from the hot spherical specimen passes through the 
aperture and is reflected by the mirror, through the IR filter, and then onto 
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the sensor head,, The emitted radiation in the wavelength range up to 1, 1 fj, 
is absorbed by the silicon cell which generates an e,m,f„ that is proportional 
to the radiation absorbed. The longer wavelengths are transmitted through 
the silicon cell so that they fall onto the underlying Lead Sulphide cell where 
they are absorbed to produce an increase in conductivity. 


Over the temperature range of interest (800°K-1700°K) the signal 

5 

dynamic range produced by the silicon cell is 10 :1 and that of the Lead 
2 

Sulphide is 10 :1„ The parameter produced by the cells that can be uniquely 
related to temperature (for a black or gray body) is the ratio of the two out- 
puts; i, e„ , 


V 


PbS 


Si 


= F (T), 


Since the ratio of Vp^g/Vg^ itself covers a wide dynamic range the signals 
from the two cells are amplified by log amplifiers and the ratio of the two 
cell outputs are fed into a differential output, so that the output signal is 
then a function of the ratio of the two outputs. Thus 


V 


PbS 

V Si 


= F (T) 


A V out = lo <S < V pbS> - lo * < V Sl > 


log F (T). 


The sensing system should be independent of the amount of energy received 
(i«e,, the size or emissivity of the hot sphere), and only dependent on the 
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amplifier. Differing time constants and thermal capacities of the PbS cell 
cell and thermistors cause an imperfect balance to exist in practice but 
the magnitude of the error introduced is not large enough to cause very serious 
inaccuracies. The cell output passes through a follower of high input impedance 
and then into a divide -by-four circuit before being amplified by the logarithmic 
amplifier. The latter has an FET, 2N3823, as the logarithmic element. 

.The silicon cell drives current into what is essentially a short circuit 

load; i.e. , the virtual ground of the 40J operational amplifier with the FET 

feedback providing the logarithmic response. The silicon cell drive current 

-10 -5 

covers the range from 10 amps to 10 amps. 

The outputs from the two logarithmic stages feed directly into a 40J 
operational amplifier operating in a linear mode. 

2. 6. 3 Circuit Setup and Tests 

The amplifying circuitry for the lead sulphide cell and silicon cell were 
tested as separate entities. Test input circuit configurations were set up and 
the logarithmic responses checked over the expected input signal dynamic 
range; the output point monitored was at the output of the logarithmic stage. 

2. 6. 3. 1 Electronic Current Tests 

The input test circuit was as shown on Figure 2.6-2. In the initial 
state the voltage at A was set to the smallest achievable positive voltage. The 
trimming resistance on the logarithmic amplifier was then adjusted so that the 


2-77 











ratio of the signals developed by the two detectors. The output of the 
amplifying electronics, in the experimental configuration, is monitored by a 
digital voltmeter, 

2, 6, 2 Circuit Description 

The complete circuitry up to the input of the telemetered or recorded 
signal is shown in zones A8-A10 on Drawing ER 47J225305, Sheet 1 in Appen- 
dix A. The lead sulphide cell is biased with a total voltage of 25,4 volts; its 
series load is not optimum in terms of signal sensing efficiency but is 
optimized for stability in the d, c, output at the cell load-cell junction. The 
load consists of two thermistors, type FA425 made by Fenwal Electronics, 
each of which has a nominal resistance of 20 KO at 22 C, In series with 
the thermistors is a 1% tolerance resistor of value 66 Kf2 and a trim 
potentiometer of 1 0 Kfi. The thermistors are bonded to the periphery of the 
PbS-Si detector combination; their function is to offset the changing resistance 
of the PbS cell as a function of temperature, A theoretical analysis of the 
circuit showed that if 


A r t R t +R l 

where ARq, £.R t are the temperature coefficients of resistance of the PbS 
cell and the thermistor respectively; R , R and R are the resistances of the 
PbS cell, thermistors and resistive load, then thermal changes around the 
sensor package should not affect the voltage applied to the input of the 
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log amplifier output was in its most stable condition with the output voltage 
at a minimum. Under these conditions the input voltage (pt A) was +1 mv 
and the output (pt B) was 0. 297 volts. 

The input voltage at A was varied from 1 mv to 1.44 volts by means of 
the variable 50 Kfi potentiometer. The relationship between input voltage 
and output voltage was plotted. The output voltage was also monitored by an 
oscilloscope to ensure that there was not an excessive amount of 60 Hz noise. 

The results of the test are shown also on Figure 2.6-2. 

The silicon cell test circuit was as shown in Figure 2.6-3. A variable 
voltage was derived from a small battery and variable resistance, and this 
voltage was fed through a 10 MO resistance into the input terminal of the 
logarithmic amplifier. The initial conditions were similar to that of the lead 
sulphide cell with the exception that the input voltage was negative. The 
initial voltages were 0.5 mv input and 0.280 volts output respectively; they 
were measured with the same DVM' s as in the lead sulphide cell test circuitry. 
In order to cover the total range of input current, the 10 MO resistor (measured 
value = 9. 945 Mf2 ) was replaced by one of 1 MQ (measured value = 1. 172 Mf2). 
The results of the test are shown on Figure 2.6-3. 

2. 6. 3. 2 Optical Checkout 

The final test consisted of irradiating the optical sensor with black body 
radiation from a 1 cm aperture of an IR Industries radiation source with a 
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Simulated Output Current from Si Cell - Amps 







set-up as sketched in Figure 2.6-4. The source was controllable in tem- 
perature over a range that extended from room temperature to about 1200°K. 


•* 8 cm 


z 


Sensor Head 






/Variable Speed Chopper 


Uj 




Industries 
Black Body 


— Black Body 

Aperture Selector 

Range of Apertures: 0.1 to 1.0 cm diameter 


Figure 2.6-4. Optical Test Set-up 


Measurements were initiated at a temperature of approximately 950°K 
and were made at 50°K intervals from this temperature up to 1200°K„ In 
addition to the output of the sensing circuitry other circuit parameters were 
monitored for each measurement made. 

The circuit was designed to make dc measurements only but because a 
simple chopper system was an integral part of the black body source, 
additional measurements were made with the input radiating chopped. The 
results are tabulated in Table 1. Plotted on Figure 2.6-5 is the output from 
the electronics with the sensor illuminated by a steady radiation. An incon- 
sistency was observed in the data for a chopped input and time limitations 
did not permit further analysis. These data were not plotted. 
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TABLE 1 


Source Temperature 


Output Volts Output Volts 

Steady Input Chopped Inpc 


974 0.157 0.14 

1026 0. 117 0. 080 

1081 0. 114 0. 070 

1131 0.065 0. 26 

1180 0. 060 0. 32 


2.6.4 Discussion of the Results 

The sensor circuit cannot be used in its present form because the change 
in impedance of the Lead Sulphide cell resulted in changing voltages being 
applied to the logarithmic amplifier that amplifies its signal output. Apparently 
the change in impedance was caused by cell temperature changes resulting 
from the proximity of the cell to the heated specimen and in spite of the inter- 
position of the IR filter. The thermistors appeared to compensate adequately 
for ambient changes. The change may have been a consequence of hysteresis 
since after irradiation the output voltage of the Lead Sulphide log amplifier 
was observed to have changed 120 mv at the highest temperature while at tne 
lower temperature (~ 1000°K) the change was of the order of 20 mv. There 
was insufficient time available to determine the exact cause of the impedance 
change. 
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The silicon cell circuitry was extremely stable, varying by only 1 mv 
over a period of about one hour. 

2. 6. 5 Re c omme ndations 

Tnere are two possible alternatives to tne present system. Tbe first is 
to use silicon cells separately filtered between 0.7 (J, and 1.1 /J. followed by a 
dc amplifying system. Tbis arrangement could be used if tne temperature 
range of interest lies above 1000°K. Tbis approach would require two 
separate silicon cells in conjunction witb a split-beam optical system. 

The second alternative is to use the current detector system of Silicon 
and Lead Sulphide but the radiation must be mechanically chopped so that an 
ac amplifying system can be used with the Lead Sulphide Cell, The electronics 
for this application would require a dc restoring system so that the logarithmic 
amplifier is not subjected to negative -going signals. 

With either of these two versions a reliable instrument can be made. 
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2. 7 SPECIMEN ROTATION MODE ANALYSIS 


2. 7 . 1 Shape Forming By Rotation of Containerless Melts 

Because of the possible importance for zero gravity shape forming by 

exploiting the equilibrium rotational forms of liquid masses acted upon only by 

surface tension and centrifugal forces, we here give a brief discussion of the 

application of electromagnetic fields for rotation mode control. The formation 

of spherodial glass, metal or other material objects by solidification while rotating 

(1 2 1 

in an equilibrium ellipsoidal form has been proposed by Deeg and others' * . 

For materials with reasonable electrical conductivity, including metals, 
semiconductors and high conductivity glasses in the molten state, controlled 
rotation can be imparted by causing the magnetic field used for position control 
to rotate at the desired rate. The levitated object will then spin up as a rotor 
of an induction motor and the rate of spin can be adjusted by suitable adjustment 
of the rate of rotation of the magnetic field. This would be implemented by 
proper phasing of the excitation to the various positioning coils. For example, 
in a six-coil system, if two opposing pairs of coils are excited at the same 
frequency but with a 90 phase shift between them, a rotating magnetic field 
will be set up at the coil excitation frequency. This type of rotational excitation 
could be applied in place of the normal excitation utilized for position control by 

Deeg, Dr. E. W., "Glass Preparation in Space," Space Processing and 
Manufacturing Meeting, MSFC, Huntsville, Ala., Oct. 21-22,1969. 

2) Frost, R. T. , "Techniques and Examples for Zero-g Melting and 

Solidification Processes" (paper presented at the Seventh Space 
Congress, Cocoa Beach, Florida, April 23, 1970), Figure 1. 
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a time sharing arrangement in which the rotational field is applied for a small 

fraction of the total time so that position control is not lost. Alternately, the 

rotational mode excitation can be superposed upon the position control coil 

excitation without interfering with the normal position control servo action. 

The latter arrangement would require a separate winding on the coupling 

transformers which supply excitation to the position control coils and would 

also require that the rotational fields not be strong as the control field. 

2. 7. 2 Dependence of Equilibrium Shape Upon Physical Parameters 

The relation between the rotational speed Oirequired to give a desired 

oblate spheroidal shape to a given mass of liquid (presumed in the molten 

state prior to solidification) can easily be calculated by considering that the 

difference in surface tension pressure due to the different curvatures at the 

pole and equator must be compensated by the centrifugal forces acting on a 

radial column of liquid in the equatorial plane. If we consider an oblate 

stable spheroidal shape with major and minor semiaxes a,b (the minor axis 

lies along the axis of rotation), the curvature at the pole can be shown to be 

Cp = -2b/a Z . The curvature at the equator can be shown to be equal to 

C e =-a/b 2 - 1/a. Assuming that fluid currents within the spheroid have been 

dissipated by viscous effects, the resulting difference in static pressures along 

a radial column of fluid in the equator plane will be given by Ap = a ( C e -C p ) 

where a is the surface tension. This pressure difference must equal the total 

accelerating force acting upon the radial column of density p (considered of 

unit cross section). The inertial f orce due to the centrifugal acceleration of 

2 

the column of fluid is given by integrating the inertial force pw r along the 
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radius, cois the angular speed. We thus obtain cr(C e “Cp) ■/ P U)2 rd 
<3^0 OJ 2 / 2. Introducing the difference in surface tension pressure between 


r = 


pole and equator, we finally obtain 


T 

a 


J \ 2 Ct 


where £ is the shape parameter b/a. The shape parameter is, of course, 

unity for a sphere and can approach smaller values limited only by stability 

considerations for rotating liquid masses. 

The equilibrium shape of a revolving liquid under its own capillary forces 

is derived, for example, in Phil. Mag., Vol. XXVIII, pg. 161, 1914. This 

derivation shows that the shape will approximate that of an oblate spheroid 

provided that the deformation is not too great. This restraint may be expressed 

2 3 

in terms of a parameter = p U) a / 8 C T. The shape will approach an exact 

prolate ellipsoid when £2 approaches unity. We see from the equation above that 

1 corresponds to that rotational speed CO which gives the shape factor £ = 1/2. 

Thus it is to be expected that the departure from spheroidal shape will begin to be 

2 3 

be apparent for such speeds. For small deformations, £= 1 - — — — — — 


For a large mass of molten material with low surface tension, we might 
consider as an extreme example the following parameters 
a = 1 0 cm 

_2 

p- 10 grams-cm 
(7 = 100 dyne-cm ^ 

4 = 1 / 2 (i. e. , minor axis =1/2 the major axis ) 
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and obtain the corresponding rotation speed 0. 3 radian per second. As another 
extreme example, we may consider a molten aluminum sphere of radius 1/2 cm 
with the same shape parameter £ - 1/2 and obtain the corresponding rotation 
speed of 157 radians per second|l500 rpm.j 
2. 7. 3 Two-Phase Induction Motor 

When the positioning coils are excited in quadrature so as to produce a 

rotating magnetic field, the resultant torque due to the magnetic dipole generated 

by the induced eddy currents can easily be calculated, provided that the skin 

depth in the conductor for the corresponding rotation frequency is large compared 

to the radius of the sphere. From the examples of required field rotation speeds 

corresponding to interesting deformations given above, it can be seen that this 

condition will normally be met for reasonable size objects. The expression for 

the torque exerted by a rotating magnetic field upon a conducting sphere 

2 3 

can be written * as T = 0. 1 (o C Irp B a per unit volume in MKS units where TJ is 

the electrical resistivity. For example, if we consider an aluminum sphere of 

-2 - 2 - 

diameter 1 cm in a field of 100 gauss ( 10 webers -m ), rotating at 100 radians/ 

-3 -3 

sec., torque per unit volume will be 4.5*10 newton-meters-meters = 0.045 

. 3 3 

dyne-cm-cm . Since the volume is approximately 1 cm , the torque will be 

about 0.045 dyne -cm. The moment of inertia of this sphere is approximately 
2 

1/7 gm cm . For these parameters we see that the ball would be accelerated 


* Zijlstra, H. , Experimental Methods In Magnetism, Yol. 2: Measurement of 
Magnetic Quantities . (Amsterdam: North-Holland Publishing Company, 
1969), p. 89. 
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at 0. 320 radians sec so that a speed of 10 sec would be reached within 

about 300 seconds. For larger spheres, the moment of inertia varies at the fifth 

power of the radius whereas the torque varies as the sixth power. Thus the 

angular acceleration in a given rotation field will increase directly as the sphere 

radius. On the other hand, the speeds required to give a given shape deformation 

vary as a ~^{^so that the acceleration time for a given deformation varies as a 

This, then , seems a practical way to achieve desire rotational shapes provided that 

the resistivity Tj is not too high. Since field strengths of thousands of gauss are 

quite practical, and since spin-up times on the order of many minutes can 

4 6 

probably be allowed, we can consider situations for resistivities 10 to 10 
times higher than that for aluminum and still achieve large deformations within 
a reasonable time. 

2.7.4 Single -Phase Induction Motor 

This method of spin-up for a solid sphere appears practical and has been 
previously studied in the laboratory. Both theory and experimentation done to 
date, however, indicate that this method is not practical for spinning molten 
spheres due to the very low torques available and the viscous torque due to 
the inevitable shape distortions of the liquid mass due to magnetostrictive forces. 
The fluid must rotate, whereas the prolate deformation of the spheroid lies 
along the direction of the main single phase field, giving appreciable damping 
torque which we estimate below. 
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2.7.5 Theory and Experiments 
2. 7. 5.1 Solid Rotor 

The torque versus speed characteristic of a single -phase induction motor 
without special starting provisions is well known and is represented schematically 
in the sketch. 

; Torque 



The torque is zero when the rotor is at rest but increases so as to accelerate 
any rotational motion, in whichever direction it is initially imparted. For 
field excitation frequencies on the order of tens of kHz which have been 
utilized for electromagnetic positioning and heating of containerless melts, the 
torque-speed characteristic can be considered linear up to the highest physically 
achievable speeds. These limiting speeds were determined in the laboratory to 
established by air drag of several thousand r pm for a 1 cm diarneter aluminum 
sphere which was simultaneously levitated and spun-up in the single phase induction 
mode. In one experiment in which air drag was eliminated by means of a vacuum 
pump, a 1 cm diameter aluminum sphere was actually burst as it reached its 
fracture stress limitation at tens of thousands of rpm. 

-2 

In the presence of air, an angular acceleration rate of 10 radians sec 

measured and appeared to be reasonably constant over the range up to 3000 
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rpm. The ultimate rotation speed was 10, 000 rpm under conditions corresponding 
to approximately 150 amperes excitation at 30 kHz into a 2 cm diameter hem- 
ispherical coil. The estimated spin-up torque was deduced as approximately 2 
dyne -cm. The near constancy of acceleration throughout the lower speed 
range indicates that the spin-up torque is not due to a torque-speed characteristic 
such as shown in the sketch above but must be due to a small quadrature component 
of the alternating field due to coil asymmetries. An auxiliary experiment was 
done in which a small flat was filed onto one portion of the spherical surface 
such that a gravitational torque on the order of several dynes-cm would be 
exerted due to the non-coincidence of the electromagnetic levitation force and 
the gravitational forces on the sphere. It was found that this asymmetry inhibited 
rotation and the sphere did not spin-up, thus giving additional confirmation to 

the order of magnitude of the spin-up torque. 

A conclusion of this work is that solid specimens should be sting mounted 
during melting under zero gravity conditions to prevent rotations where these 
are not desired. Spin-up to bursting speed of unmounted specimens appears to 
be inhibited by air drag and this inhibition would remain provided inert gas 
environments are used at reasonable pressures. 

2. 7. 5. 2 Liquid Rotor 

In a number of laboratory levitation experiments in which aluminum 
spheres of 1 cm diameter were melted and superheated, no rotation could be 
observed. This is attributed to the rather large deformations from spherical 
shape due to the combination of magnetostrictive and gravitational forces. It 
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remains to examine whether significant rotational motions can be imparted as 
a result of single-phase induction motor action when molten spheres are being 
processed by electromagnetic fields in the weightless condition. 

In a zero gravity environment, both the spin-up torque due to single -phase 
induction motor effects and the asymmetry and consequent viscous damping due 
to magnetostrictive forces will be proportional to the exciting field strength. 

2. 7. 6 Shape Deformation Due to Magnetostrictive Forces 

The asymmetry induced in a conducting liquid sphere due to magneto- 
strictive forces can be estimated in the following way. Computations of elect- 
romagnetic force fields within conducting spheres located in uniform oscillating 
electromagnetic fields were summarized in the reference * . These show 

that the magnetostrictive forces act radially inward around the equatorial 
regions of the sphere and vanish toward the axis of the sphere (assumed in the 
direction of the magnetic field) as well as diminishing at higher latitudes. 

Such a force field will result in a circulation of fluid which is radially inward 
in the equatorial regions and poleward along the axis (see sketch). The 
resulting flow system consists of two ring vortices between which there will 


0 

C 



* Frost, R. T., et al, Field Management for Positioning and Processing of 
Free Suspended Liquid Materials. General Electric Company, Contract 
NAS8-24683, Modification No. 2, Task IV, Final Report, May 15, 1970. 
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be a repulsive effect, leading to deformation of the sphere into a prolate shape 
elongated in the magnetic field direction. Although the fluid dynamics problem 
has not been solved taking account of viscosity and surface tension, an order 
of magnitude estimate for the deformation can be obtained by applying Bernoulli' s 
theorem between the points a, b and c. Since these three points are stagnation 
points, the corresponding kinetic energy per unit volume will vanish and we can 
write p - p = W, where W is the work done by the electromagnetic force in 
moving a unit volume of fluid radially inward from a to b against the pressure 
gradient due to the nonuniform curvature of the deformed surface film. Part 
of the electromagnetic work will go towards the overcoming of viscous forces 
but this will be neglected in this order of magnitude estimate. In this 
approximation, the pressure at points b and c will be the same because both are 
stagnation points and there are no electromagnetic forces acting along the 
axis of the spheroid. The pressure differential can thus be taken as due to the 
difference in surface tension pressure at the points a and c as derived earlier 
in this section to obtain 

P b " Pa = ( 1 + ~J2 ~ 2 

The work done by the electromagnetic forces per unit volume of fluid moving 

from a to b can be written as 

a 

r > > > 

W = - / ( j x B) . dr 
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where j is the current density and B is the magnetic flux. Both j and B will 

die off rapidly as we proceed inward from point a and these variations were 

plotted in the reference $ For our present order of magnitude purposes the 

effective B can be taken to be the magnetic field at the spheroid surface and 
a 

the integral f jdr can be equated to the total eddy current flowing around 

Jo 

the equatorial region of the sphere per unit latitudinal distance. 

W = B / jdr = l u B 

The total circulating eddy current induced in the spheroid can be found 
by equating the magnetic dipole moment derived by Smythe ** to the 
classical expression for the dipole moment in terms of current times area. 

This results in the identity 

3 2 

71 a G (x)H = 7T a i 

or 

i = aG(x)H 

For a representative laboratory situation in which a 1 cm diameter aluminum 

sphere is levitated and melted in a field of the order of 100 gauss @ 30 kHz 

the resulting current is approximately 50 amperes and is distributed over 

approximately 1 cm of latitudinal distance. G(x) = 1. In gaussian units the 

-3 

total work W = i. B = 5 emu current 100 gauss = 500 ergs cm . Equating 

this to the difference in surface tension pressures at point a and point c gives 

for the deformation parameter £ the relation 

_!L ( 1 + — - 2 £ ) = — F ( £ ) = 500 ergs - cm 

a ^ 2 a 


Frost, R. T. , et al , Field Management for Positioning and Pro cessing of Free 
Suspended Liquid Mate rials. General Electric Company, Contract NAS8-24683, 
Modification No. 2, Task IV, Final Report, May 15, 1970 

Smythe, W. R. , Static and Dynamic Electricity (New York: McGraw-Hill 1950), 
p. 398. ” 
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giving a deformation parameter of about 0. 95. 


2.7.7 Viscous Damping Torque 

The single -phase induction motor torques are perpendicular to the magnetic 

field direction and hence will tend to spin the spheroid up in a direction normal to 

its axis of elongation. Since the direction of elongation remains fixed in the 

direction of the alternating magnetic field, viscous damping effects will occur 

when rotational motion is imparted to the entire fluid mass as a result of single 

phase induction to torques. We can estimate the limiting speed for rotation 

set up in manner by equating the spin-up torque derived for the solid sphere 

experiments described earlier to the opposing torque due to the viscous 

dissipation when the fluid mass rotates but maintains its elongation in a fixed 

direction. 2 n T= E, where E is the viscous energy dissipated in a single 

rotation. The viscous energy dissipated in 1 / 4 rotation can be estimated as 

approximately twice that due to deformation of this sphere from the prolate 

into a spherical shape. The viscous decay problem has been studied for the 

latter case in terms of a viscous damped sphere initially defoormed in a 

prolate fashion * . For this model, the viscous decay time can be written 
2 

as t^- ywhere V is the kinematic viscosity. Estimating that the change in 

surface area of the spheroid deformed with a shape parameter 0. 95 is of the 

2 

order of 0. 0015 cm and taking a surface tension of 900 dynes per cm, the 
additional surface energy of 0. 0015*900 = 1.4 ergs would be dissipated in a time t 

* Lamb, Sir Horace, Hydrodynamics (New York: Dover Publications, 1945), p. 640. 
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of the order of one second for spheres of many liquid metals of the size con- 
sidered here. * Considering that this energy would be dissipated in approximately 
1/8 of a total rotation for the analogous problem of the rotating fluid mass con- 
stantly deformed in a given direction, the power dissipation due to viscous forces 
will be 8 x 1.4 ergs/ 8 seconds corresponding to a rotational period of 8 seconds. 
Since the viscous power dissipation will increase as the square of the angular 
velocity to we can write for the rate of viscous dissipation T CO the product 

v 2 

T co = 1 lo Z er & S x | =2.6 ergs-sec -a. , where CO is the angular 

8 sec yco^ J ° 

speed corresponding to a rotational period of 8 seconds. Finally, since the 

laboratory experiments referred to indicated spin-up torques on the order of 

2 dyne-cm, we conclude that the maximum rotational speed which can be imparted 

to the viscous, deformed spheroid will be approximately 1 radian per second. 

At lower field strengths such as in the case where the spheroid would be allowed 

to solidify under the influence of relatively smaller field strengths, both the 

deformation and spin-up torque would fall off proportionately so that the same 

argument would apply. Our order of magnitude estimate for the very low angular 

rates which can be induced by single phase induction action in suspended liquid 

metal systems is confirmed by the failure to observe any appreciable rotational 

motion for molten levitated aluminum spheroids in the laboratory. 


* Frost, R. T. , "Techniques and Examples for Zero-g Melting and Solidification 
processes, " (paper presented at the Seventh Space Congress, Cocoa Beach, 
Florida, April 23, 1970), Fig. 1. 
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3. 0 MODIFICATION OF SIX-COIL DROP TEST UNIT 

3. 1 MODIFICATION TO INCORPORATE TECHNICAL DEVELOPMENTS 
When the drop test unit was received it was tested and found in good 
working order, therefore no repairs were necessary. Figure 3. 1-1 is an overall 
photograph of the equipment as modified. The modifications made as a result 
of technical developments of EMLS were as follows: 

« Power Amplifiers (Figure 3.1-2) 

The operating frequencies of all six channels were changed to those 
near 100 KHz. The multipliers were modified to go from the 
original heating mode (all amplifiers delivering maximum power 
with no specimen position error) to one where the amplifiers 
only deliver power for positioning. 

& Pre -Amplifiers 

Six pre -amplifiers were added to utilize the position sensors 
developed under this contract. These pre -amplifie rs were mounted 
underneath the coil assembly base within a shielded housing, 
o Position Sensors 

New position sensor assemblies were built for the drop package 
using the sensors and approximate optical system (but without 
mirrors) of the prototype system described in the Sections above. 

• Servo Electronics (Figure 3. 1-3) 

Modification in gains and break frequencies were made to be 
compatible with the new coils, frequencies and specimens. 
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igure 3. 1-3 Differential and Servo Amplifiers Position Sensing 





» Coil Assembly (Figure 3. 1 -4 ) 

New coils and tank circuits and a new lighting system was built on the 
original coil assembly base as shown in Figure 3. 1-4). Also added 
were lights, (on the left corner of the assembly pictured in Figure 
3. 1-4) six of which indicate which power amplifiers is driving and 
a seventh which indicates an overtemperature kick-out of one or more 
power amplifiers. A small regulator for the specimen lighting system 
was added to stabilize the unregulated voltage and maintain a quiet 
light output. Some interaction with the servo had been found as the 
light intensity changed according to the servo drive requirement. The 
regulator eliminated this interaction. The prototype system described 
in this report defined a sophisticated infrared signal positioning 
scheme but the drop package uses a painted white ball, the entire coil 
enclosure being flat back. 

3. 2 ADDITIONAL CHANGES 

The modifications in addition to the technical developments of EMLS were 
as follows: 

• Power Amplifiers 

Over temperature sensors (thermistors) were mounted on each power 
transistor assembly. A board of electronics (seen on the right hand 
end of Figure 3. 1-5), containing control circuitry and relays in the 
power lines of each power amplifier, was located on the power 
amplifier assembly. This board also has inputs from thermistors 
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Figure 3. 1-4 Positioning Coil Assembly 



Figure 3. 1-5 Power Amplifier Assembly with Overtemperature 

Protective Circuits 
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mounted on each coil and so provides protection against over tem- 
perature in either the coils or power amplifiers. 

A small board containing three additional oscillators was mounted on 
the power amplifiers assembly and can be seen on the end of the 
assembly on the right side of Figure 3. 1-1 • 

© Servo Electronics 

A limiter was placed on the position error line, but not the rate 
line, into the servo amplifier. The purpose of this limiter is to 
prevent large position errors from driving the power amplifiers at 
full power. They can be readily removed if full force as a result of 
a position error is desired. Because full force is available for damping, 
the time to show convergence from a position saturation condition is 
shortened. This is important for demonstrating damping in the MSFC 
drop tower facility where the specimen rests directly over the 
bottom coil at the instant the drop starts. Without position error 
limiting this specimen will be accelerated by a force of about 200 
dynes for the first centimeter or so, and will overshoot the center 
with a rather high kinetic energy. The time to damp out this large 
energy is much longer than that required with position error limiting 
which limits the coil force to about 50 dynes. Note that with position 
error limiting there is no rate error limiting so damping has the full 
200 dynes available near any coil. 
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There has been no change in mounting centers, outline dimensions, 
or cable lengths in the refurbished drop test unit. 

Two copies of a technical manual will be delivered with the drop test unit. 
This document will be the same type as was sent with the original equipment* 

It will include the following items: 

• Description of operation utilizing a block diagram 

• Operational and adjustment procedures 

• Circuit diagrams (schematics) 

• Photographs of circuit board layouts and construction 

• Position sensing error curves 

• Telemetry calibration data 

• Wiring and cabling diagrams. 
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APPENDIX A 


ELECTRONIC BLOCK AND SCHEMATIC DIAGRAMS 

Included in this Appendix is the Block Diagram to implement the functions 
of the Electromagnetic Levitation System (EMLS), General Electric Drawing 
No„ ER 47D225301. The detailed electronic circuits of the blocks enclosed by 
the broken line in the upper left hand portion of this drawing are depicted on 
Sheet 1 of General Electric Drawing No. ER 47J225305 which follows the Block 
Diagram. It should be noted that these circuits represent the first effort to 
mechanize the design equations for determining drive currents as a function 
of position errors generated by the position detector and servo system. These 
circuits include the elements necessary to sense position from information 
supplied by the IR sensor to accomplish the coordinate conversion, to provide 
servo lead compensation and to compute the driving currents required. The 
circuits as shown cannot be considered to be fully designed. No breadboards 
have been made to permit testing in the laboratory nor have the required analyses 
been performed. Therefore, these circuits may not represent those which would 
have been utilized had not the effort been terminated before the design was 
complete . 

The circuits to mechanize the blocks shown on the right hand side of 
drawing ER 47D225301 are detailed on Sheet 2 of Drawing ER 47J225305. No 
discussion of these circuits is included in the present report but they were 
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reviewed in detail in Section 2. 4 of the Final Report dated June 15, 1971 for 
Contract No. NAS 8-26157 entitled, "Free Suspension Processing Systems for 
Space Manufacturing. " 

The blocks shown below the broken line on Drawing ER 47D225301 outline 
the conceptual mechanization for a system to capture the spherical specimen 
after melting and re solidification was completed. No mechanization of this 
capture system had been effected at the time of the contract termination and no 
schematic diagrams are presented. 
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Explanatory Notes on Program OPS 


Lines 00070 
« 

100 

HO 

150 

170 

260 


LLi, statement functions for the four unit vectors originating at 
(0, 0, 0) and directed toward the center of the coils. 

FFi, statement functions for the solution of the unit force 
vectors. 

K = 17. Spring constant - set to allow maximum current at 
point in field. 

Output description 


360 


400 Position of point specimen on end of string. 

i 

I 

420 

430 IVEL, VMAG: pick the velocity unit vector off string and the 

j magnitude . 

5&0 


570 Start up and solution for first set of forces desired 

i 

1200 

1210 Find the current at each coil to deliver the desired forces; find 

13’40 the delivered force components. 

1370 Let specimen move according to delivered forces and continue. 

I 

l 

• 

2000 
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\ HEAD'/ 

‘ SLI ST OPS* CURENT* CNVRTX* TETMFO* MAGMFO* CMVRTF 
05/ 1.5/72 09.000 


00010 * OPS PROGRAM MFO 12/28/71 

09020 * 

03030 REAL LL1 *LL.2*LL3*LL4*L1*L2*L3*L4*L1 1 * L2 1 * L3 1 * LA 1 * 
000A0 & MAGP* MAGV* MAGF* 1< 

03050 & 11*12*13*1 4* MAGFF 

33060 * L IS DIRECTED TOWARD COILS 

00070 LL1CX* /*Z> = +0.816492*/ -0.57735*2 

03080 LL2CX*/*Z)=+0.816492*X +0# 57735*2 

00O90 L.L3C X* /, Z > = "0.8 1 6492*X +0*57735*2 

03130 LL4 ( X * / * Z ) = -0. 816/192 + / -0.57735*2 

001113 * F IS DIRECTED TOWARD ORIGIN 


00120 

FF1CX* y*z> = 

+0. 612373*/ 

-0. 433013*2+. 25 

00130 

FF2CX* /*Z)« + 3. 6 1 2 3 7 3 * X 


+ Q. 43 30 1 3M-Z + . 25 

03140 

FF3CX* /* Z)»-0. 612373+X 


>0. 433013*2+. 25 

03153 

F F4 C X * / * Z > = 

-0.612373*/ 

-0.43301 3+Z-t .25 , 

00163 

S S ( A * B * C ) = S Q R T < A * * 2 + B * * 2 + C * * 2 ) 

I 

00173 

K= 1 7 • 


i 

03130 

C 1=0.8 1 6492 


1 

00190 

C2=0.. 57735 


i 

00200 

C7= 1 . 33 



002 13 

C6=SQRKK*C7 + 4. > 


5 

002 23 

C8 = 0 * 2 




i 032 30 * C 6 DAMPING CONSTANT* C7 MASSCGR)* C8 PRINT INTERVAL ! 

i 002/0 * RESERVED FOR VOLTS TO X*V*Z» CONVERSION ■'! 

! 03250 * l 

j 03260 PRINT 888 + j 

• 00270 838 FORMATC ( 1 5X* ‘’OUTPUT FORMAT" ) / i 

: 00280 & ( 1 5X* M DC“X* /* Z*. COMPONENTS * MAGNI TUDE" > / . 

: 00290 ' & ,( 1 5X* " M=MAGNITUDES: CM OR D/NES OR AMPS")/ j 

00300 &// 

| 03310 & (2X*"TIM£"*9X*"-1*SPECI MEN POSITION DC")/ 

! 00323 & < 13X*"SPECIMEN VELOCIT/ DC")/ 

! 00333 & ( 2 IX* "DESIRED FORCES DC")/ i 

! . 00340 & ( 17X* "DESIRED COIL FORCES M">/ ; 

: 00353 & ( 23X* "CURRENT INPUT M")/ 

1 00360 : & < 1 9X* "DELI VERED FORCES DC")///) j 
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00370 
00380 
00390 
00/100 
00410 
00420 
2043 0 
004 40 
004 50 

004 60 
034 7 0 

0048 0 

0049 0 
00500 
00510 
00520 
03522 
00 524 
0052 6 
00533 
00540 
00 550 
00560 

05 570 
00530 
03 590 
03600 
03610 
33620 
03630 
00640 
03650 
00660 
00670 
0068 0 
00690 
007 00 
007 10 
00720 
00730 
00740 

007 50 
00760 
00770 
0078 0 
00790 
038 00 

008 1 0 
00820 
83o 3 0 
008 40 
008 50 
008 60 



'I 6 


~>7-‘ 


I VEL 

TO COIL. 3 CNTR 
VZ=-0. 577353 




LIME ( I VEL = 1 ) 


PRINT 35; PRINT 35 

* POSITION OF STING 

PRINT 35 
X=-0. 67400 
Y=-0. 51845 
Z= 0.8 4321 
READS IV EL • VM A 6 
GO TO 91*92*93*94)* 

* V UNIT VECT PARALLEL 

91 VX = 0 • 3 1 64 9 6 IVY=0. 

GO TO 99 

* v UNIT VECT FROM STING TOWARD ORIGIN CIVEL=2> 

92 VX = 0. 56233 JVY=0. 43294 ; VZ=-0. 70414 
GO TO 99 , 

* V UNIT VECT NORMAL TO STING CIVEL = 3>' 

93 VX-0. 09123 ;Vf=0. 89397 *VZ=-0. 43867 
GO TO 99 

* V UNIT VECT* -OFF STING DUE TO VOUCH FORCES 

94 VX= 0 • 8 432 ; Vr'--0.4672 ; VZ=-0«2658 
99 VX= VX*VMAG ) VY=VY*V MAG J VZ=VZ*VMAG 

* INITIAL CONDITIONS: SPECIMEN AT X*Y*Z WITH VELOCITY 

* VX> V Y* VZ AFTER DT SECONDS IS AT X1*Y1*Z1 WITH 

* VX* VY* VZ: COIL ACTION BEGINS MOW 

T=0.0 
DT- .01 

MAGP-SSCX* Y*Z> 

MAG V- SSC VXi VY* VZ ) 

IFCMAGV.EQ.0.0) GO TO 102 

* V UNIT VECTOR COMPONENTS 

VXU= VX/M AGV 
VYU= VY/MAG V 
VZU= VZ/MAGV 
' GO TO 4 

VXU= 0 . i VYU= 0 . ; VZU= 0 . 

PRINT 5 

FORMAT < IX* "INITIAL CONDITIONS 
PRINT 3 0 * X * Y * Z * M A G P 
PRINT 30* VX* VY* VZ*MAGV 
PRINT 6 


ci 7 '■'A 


i 

! 

\ 
? 

VELOCITY 


102 

4 

5 


T=0.00”> //) 


DISPLACEMENT 


6 FORMATC 1 X* 

♦ COMPUTE INITIAL 

X 1 =X+ VX*DT 
Yl= Y+VY*DT 
Z 1=Z+VZ*DT 
MAGP=SS(X1* Y 1 * Z 1 > 

I F( MAGP. EQ. 0. 0) GO 

* -P UNIT VECTOR 

X1U=-X1/MAGP 
Y 1 U=- Y 1 /MAGP 
Z 1 U=-Z 1 /MAGP 
GO TO 10 

103 xiu=0 .; yiu= 0. ; z iu=0. 
♦END OF START UP SEQUENCE 


//) 


TO 1 03 
COMPONENTS 
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00870 
00880 
00890 
00900 
009 10 
03920 
00930 
00940 
03950 
, 009 60 
00970 
00980 
03990 
01000 
01010 
01023 
01030 
01040 
01050 
01060 
01070 
01080 
01390 
01100 
01110 
01120 
01130 
01 140 
01 1 50 
01160 
01170 
01 180 
01 1 90 
01200 
01210 
01220 
01230 
01240 
01250 
01260 
01270 
01280 
01290 
01300 
01310 
01320 
01330 
01340 
01350 
013 60. 


1 0 CONTI NUE ' ’ " — • 

T= T+DT 

* COMPUTE L’SCP PROJECTIONS ON COIL VECTORS) 

L 1 =LL 1 ( Xj Y* Z ) 

L2 = LL2( X.» Y> Z) 

L3=LL3(X> Y*Z) 

L4=LL4( X, ■{> Z> 

LI 1=LL1CX1* Y 1 * Z 1 ) 

L21=LL2CX1>Y1>Z1) 

L31=LL3(X1> Y1>Z1 ) 

L41=LL4<X1*Y1>Z1) 

* COMPUTE F* S 

F1 = FF1(X1> Yl> Z1 ) 

F2=FF2(X1 > Yl> Z 1 ) 

F3 = FF3(X1.» Z1 ) 

F4=FF4( X1*Y1»Z1) 

* LOOK FOR SMALLEST ~F 

ADD 1 = AM I M 1 C F 1 < F2 > F3 > F4 ) 

* COMPONENT VELOCITIES ALONG COIL AXES 

DL 1 DT- <L1 1 — L 1 )/ DT 
DL2DT- ( L2 1 -L2 ) /DT 
DL3DT=(L31-L3)/DT \ 

DL4DT = ( L4 1 -L4 ) / DT 

* ADJUST F’S 

F 1 - K* ( F1-ADD1 ) + C 6*DL 1 DT 
F2=K*( F2-ADD1 ) +0 6*DL2DT 
F3=K*( F3-ADD1 )+C6*DL3DT 
F4= K* C F4-ADD1 ) + C6*DL4L')T 

* CHECK AND ADJUST FOR -F’S 

ADD2 = AMIN l(Fi*F2*F3*F4) 

F1=F1 -ADD2 
F2= F2- ADD2 
F3=F3-ADD2 
F4= F4- ADD2 

* g.\T D OF F’S CALCULATIONS 

* CALL SUBROUTINE FOR COIL CURRENTS-”CUnENT” 

CALL CURENTC FI > F2, F3> F4» L 1 1 » L2 1 >L3 1 > L4 1 > I 1 > 1 2 > I 3> 1 4 ) 

* CALL SU3RT. FOR CONVERSION OF X1«Y1»Z1« TO 

* X1P.» Y1P/Z1P, -GRAY’ S ORIENTATION OF THE TETRAH. 

* "CNVRTX" IN WHICH GRAY’S COILS A,3,C>D ARE OPS’ 

* COILS 4 > 3 > 1*2 RESPECTIVELY 

CALL CNVRTX ( X 1» Y1>ZLX1P< Y 1 P> Z 1 P) 

* CALL SUBRT • "TETMFO’’-MODI FI F.D "TETRA OF GRAY- - 

* TAKES X 1 P> Y 1 P* Z 1 P AND CURRENTS AND COMPUTES FX* FY* FZ 

* IN GRAY'S COORDINATE SYSTEM 

CALL TETMFOCXl P, Y 1 P> Z 1 P» I 1 * I 2> 1 3> I H, FXP* FYP> FZP) 

* CALL SU3RT • "CNVRTF”- CONVERTS GRAY’S FORCE COMPONENTS 
*<FXP».».) TO OPS COMPONENTSC UNIT VECT. FORM) FFXU> . . . 

CALL CNVRTFC FXP, FYP> FZP* FFXU* FFYU> FFZU* MAGFF) 
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r 01 3 70 * COMPUTE DESIRED FORCE VECTOR COMPONENTS 


01380 


FX=0.0*F1 -C 1*E2 +C1*F3 +0.0*F4 

J 

| 


01390 


FY=-C1*F1 +0.0*F2 +0 * 3i‘F3 +C 1 *F4 

{ 


01/100 


FZ=+C2*F1 ~C2*F2 -C2+P3 +C2*F4 

j 


1 01413 


MAGF=SS< FX, FT, FZ ) 

i 

\ 


; 01420 


IFCMAGF.EO. 0.0) GO TO 104 

4 


' 01430 

* F 

UNIT VECTOR COMPONENTS 

a 


: 01440 


FXU=FX/MAGF 

j 


01453 


FYU-FY/MAGF 

') 


1 01460 


FZU=FZ/MAGF 

" 1 


; 01470 


GO TO 7 

i 


01430 

1 0 4 fx u= 0 . ; f v u- 0 • ; fz u- 0 . 

4 

i 


; 01491 

7 

IFCT.GT.0.0) C8= .01 

$ 

■1 

! 


. 0 i 4 9 2 


I F( T . G T . 0 . 02 ) C8» • 02 . 



: 01494 


I F C T » G T • 0 * 1 > C3-0.2 i 



: 01496 


IF<T«GT» 1 • 5 ) C 8 = 0.5 ! 



■ 01498 


TM=T/C8 • 



. 01500 


TRT= A I N 1 ( TM + .035) -j 



[ 01510 


TEST=ABS< TRT-TM) 



01 520 


I FC TEST. LI:V. 0.0001 > GO TO 20 


•' ? 

•> 

: 01 530 


GO TO 15 ; 


i 

i 

01540 

20 

PRINT 25, T,Xl: , Ti; , Z T , t»1AGP 

1 


01550 

25 

FORMAT ( 2X, //, F8 .2, 1 P4E14. 3 > 


i 

. 01560 


PRINT 33, VX , VY , VZ:_ , M A G V 


j 

i .01 570 


PRINT 3 0 , FX > FY. , FZ. * MAG F 


* 

: 01580 


PRINT 3S,F1,F2,F3,F4 


i 

; 01590 


PRINT 30, l 1 .* 12 . T 3, T 4 


M 

1 

01600 


PRINT 30, FFX , FFY , FFZ MAGFF 


; 01610 

30 

FORMAT < 5X, 1 P4E14* 3) 


1 

j 

! 01620 


PRINT 35 

' 

, j 

; 01630 

35 

FORMATC / ) 



; 01640 

* 




: 01650 

15 

CONTINUE 


.i 

: 01660 

*CO. 

•IPUTE ACCELERATION COMPONENTS DUE TO 

ACTING FORCES! 

01670 


AX=< FFXU+MAGFF) /C7 


« 

: 01680 


A Y= < FFYU*M A GFF ) /C 7 



• 01693 


AZ= C FFZU*MAGFF)/C7 



: 01700 

* TRANSFER OLD. POSITION AND COMPUTE NEW 



j 01710 

* 

TRANSFER 



j 01720 


X=X1 


* ' 

i 01730 


Y= Y 1 



; 01740 


Z=Z1 


j 01750 

* 

NEW POSITION 



01760 


X1=X1 +( VX + 0. 5*AX*DT> *I)T 



1 01770 


Yl=Yl+( VY+0. 5*AY*DT>*DT 



[ 01780 


Z 1 = Z1 +(. VZ + 0. 5+AZ *DT>*DT 



i 01790 


MAGP=SS(X1, Yi,Zl) 



[ 01800 

„ - 

I F< MAGP. EQ. 0.0) GO TO 106 



*.4 i 
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r 01810 
! 018 20 
; 01830 
: 018 40 
; 018 50 
018 60 
: 0 1870 
i 01380 
: 01890 
i Cl 903 

; 01910 

i. 01920 
i 01930 
L 01940 
i 019 50 
| 01960 
| 01970 
!■ 01930 
' 01990 
02000 


’" 00010 

! 00020 
; 03030 
; 00340 
: 03S50 
: 00060 
: 03070 
; 00080 
! .00390 
00100 
j 00110 
| 00120 
' 00130 
: 00140 

• 00150 
| 00 1 60 
[ 00170 
{ 00180 
I 00190 

l 00200- 
! 00210 
i 00215 
00220 
; 00230 

j 00240 
j 00250 
; 00260 
: 00270 
! 00280 
: 00290 
. 03330 

; 00310 

• 00320 


P UNI T VEC TOR “COMPONENTS' FOR -P 
X1U=-X1 AM AGP 

yiu=-yi/magp 
Z1U=-Z1/MAGP 
GO TO 8 

106 xiu=0.; yiu=0«;ziu=0» 

* NEW VELOCITY 
8 VX= VX+AX+DT 

VY=VY+AY*DT 
VZ=VZ+AZ*DT 
k A G V=SSC VX > VY t VZ ) 

I FC M AG V • EQ • 0 • 0) GO TO 103 

* UNIT VECTOR COMPONENTS FOR V 

VXU= VX/MAGV 
VYU= VY/MAG V 
VZ U= VZ /MAG V 
GO TO 10 

103 VXU=0»i VYU=0. I VZU=0. 

GO TO 10 

END ■ ' ■ . 



• J 


i 

■» 

j 

, \ 
i 
i 

X 

l 

■4 

I 

l 

4 

J 

j 


1 

•-» 

1 


" F 1 > F2 , F3 , F4 , XL 1 1 , XL2 1 , XL 3 1 > XL4 1 >' 

& XI \> XI2>XI3>XI4) 

XI ICA,B,O = C0.5*ALOGCA)-1 .2464*B + C> ! 

i 

1 

I FC FI .GT.0.3) GO TO 10 i 

XIl=l.E-6 ’ 

GO TO 12 i 

10 X I 1 = EXPC XI I ( FI * XL 1 1 > 2« 7764 ) ) j 

12 IF(F2.GT.0.0> GO TO 20 ! 

XI 2= 1 • E~ 6 | 

-GO TO 22 . l 

20 X I 2= EX PC XI I C F2* XL 2 1 .» 2« 7733) > I 

22 IFC F3. GT.0.0) GO TO 30 j 

XI3=l.E-6 j 

GO TO 32 ; 

30 XI 3=EXPC XI IC F3*XL3 1*2. 761 1 > > ] 

32 IFC F4. GT.0.0) GO TO 40 I 

X I 4 - 1 « E- 6 j 

GO TO 42 | 

40 XI 4= EXPC XI I < F4* XL4 1*2* 76 63 >) ] 

42 CONTINUE j 

* ADJUST CURRENTS J DO NOT EXCEED MAX BUT KEEP RATIO 

XMAX= 45 • 1 

XTRY=AMAX1 CXI 1 > XI 2* XI 3> XI 4> ' j 

I FCXTRY.LE.XMAX) GO TO 50 j 

FACT=XMAX/XTRY I 

XII =XI 1 *FACT ’ j 

XI 2=XI 2* FACT ] 

XI 3 = X I 3 * FAC T j 

XI 4 = X I 4* FACT I 

50 CONTINUE J 

RETURN A 

END ‘ ; 



r 00010 
; 00020 
00030 
> 03040 

; 00050 

| 00060 
j. 03870 
03380 

*" SUBROUTI NE C NVRTX ( X 1. Y 1 * Z 1 i X I P. Y 1 PV Z IP) - 

A= S wiiT ( 1 * / 3 * ) 

B= SQRTC 2 • /3 . ) 

X1P= A*X1 + 0 » 0* Y 1 + B*Z 1 

Y1P=0.0*X1 +1.0*Y1 +0.0+Z1 

Z1P= -B*X1 +0.0*71 + A*Z 1 

RETURN 
END 

800 10 

SUBROUTINE TF 

;TMFO( X. Y. Z. Z 1. Z2. Z3. Z4. SUMX. SUMY. SUMZ) 

I , 03823 

DIMENSION EYAC 4) . E73C 4) . EYCC 4) .EYDC 4) . SKC 4 ) . FXC 4 ) . FYC 4) . FZC4 i 

1 00024 

I FC X • EQ . 0 • 0 ) 

X= 1 . E- 6 

i 00026 

IFC Y.EQ.0.0) 

Y= 1 . E- 6 

| 00028 

I FC Z * EQ .0.0) 

Z=1 .E-6 

; ■ 00070 

A= 1 . 1 5 


| 00080 

DO 10 1=1.4 


i 00090 

EYAC I ) = 1 * E" 6 


f, 00100 

EYBC I ) = 1 • E” 6 


; 00110 

EYCC I ) = 1 • E“ 6 


00120 

EYDC I ) = 1 • E- 6 

1 • 

00130 

10 CONTI NOE 


00 1 40 

EYCC 3 ) = Z 1 


88150 

ETDC 1 ) =Z2 


00160 

EY3C 4 ) = Z3 


30170 . 

EYAC 2 ) = Z4 


j 00 1 8 0 

RS=0 » 5 


; 80135 

SCM=3 .05 


j 03190 

S=SCM*.0i;X= 

X*.0i; Y=Y+.0UZ=Z+.01 

; 08208 

AM=A+ . 3 1 



00210 SK( 1 ) = RS/0« 0477 

00220 SKC2)=RS/0.0436 

30230 SKC3)=RS/0.0494 

00243 SK(4)=RS/0.0406 

00250 RS=RS*» 0 1 

00260 PI=3. 1415926 

00270 F1=0. 5*SQRTC 3* 0) 

00280 G 1 = C 1 . 0/6. 0) *SQRTC 3.0) 

00290 F2=3. 5 

00300' F3= 1 • 0/6* 0 


■ 003 1 0 

| 00320 

j 03330 
i 00340 
! 30350 

t 00360 
1 00373 

■ 03383 
00390 
03430 

; 034 1 0 
j 03420 
j 00430 


F4=C 2. 0/3. 0)*SQRTC 2. 0) 

G2= C 1 . 0/12.0)*SQRTC 6.0) 

F5= C 1 . 0/3.0 )*SQRTC 2.0) 

F6= C 1 .0/3. 0) +SQRTC 6. 0) 

F7=l. 0/3*0 
DO 59 J=l/4 

XA=F1*CX+G1*S)-F2*C Y+F2*S) 

YA=F3*C X+G1*3)+G1*C Y+F2*S) -F4*C Z+G2*S) 
ZA=F5*CX+G1*S)+F6*C Y+F2*S ) + F7* C Z + G2* S ) 
ZAC=ZA 

RH 0 A = S Q R T C X A ■ * X A + Y A * Y A ) 

PHI A=ATAN2C YA.XA) 

BPHI A=0 . 0. , , . . T.. ..- P.'.-. , 
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r 00440-' 
00450 
!' 00463 
j 03470 
’ 0343 0 

: 0349 O 

j 00503 
00 5 i 0 
| 00590 

| 03530 

■ 00540 

| 00550 

I 03560 
i 03570 
| 03580 

j 33590 
00600 
| 03610 

| 0069.0 

l 03638 
| 00640 

i 03650 
036 60 
; 00670 

; 30630 

j C0693 

? A r> r% z’* 

;• 'OOtVU 

j 00710 
j 03793 
03730 
| 0G740 

} 03750 

j 00760 
j G0770 
j 0078 0 
I 00790 
| 003 00 

| 003 1 0 

i .036 2 0 
| 00330 

{ 00840 

| 003 50 

| 003 60 

j. 003 73 
| 03880 

! 00890 

j. 009 00 
| 00910 

; 00990 

i 00930 
!■ 00940 

i 00950 
009 60 
! 00970 

i 03980 



CALL MAG ( J* AM/ EYA, RHOA/ ZAC/BRHO/ BZAC > 

BRA1=BRH0 

IFCRHOA.LT. 1.0E-05) BRAT =0.0 
BZ A 1 = BZAC 


SIH=(EXP( 2.S*SKC J) ) - KXPC - SKC J) ) ) /2 • 0 J 

SIK1=C EXPC SK( J) >-EXP( “SKC J) ) 1/9.0 ! 

Q-- C SIH-SIMC 2* 3 + SKC J) ) ) /( SI HI *SIH 1-+SIMC SKC J) > *SINC SKC J) ) ) 
GX= 1 . 0-C 0. 75/SKC J ) ) *Q j 

DRMOA* • 000 ! j 

Erl 0 A= RH 0 A + DRH 0 A s ] 

CALL MAGC J/ AM/ ETA/ RHOA/ ZAC/BRHO/ BZAC) j 

BRA8=BRH0 j 

BXA2=BZAC ■ i 

DBURHO= C BRAS -BRA 1 ) /DRHOA \ 

RH 0 A = RH 0 A - D RH 0 A ' ] 

ZAC=ZAC+. 0001 j 

CALL MAGC J/ AM/ ETA/RHOA/ZAC/BRHO/BZAC) j 

BRA3=BRH0 j 

BZ A 3= BZAC ' " -j 

EM U= 4 . 0 * P I * 1 . 8 E- 0 7 j 

D3DZ= C 3ZA3-BZA1 ) / 0 • 030 1 ] 

D3RDZ=(BRA3-BRA1 )/0. 0031 | 

DB D Z R= < B Z A 2 - B Z A 1 > / 3 . O 301 1 

BT= 5QRTC BRA 1 *BRA 1 +BZ A 1 *BZA 1 ) j 

SA=BZA1/BT ] 

C A~BRA 1 /DT j 

DB T D S : = S A * C A * DB DZR+SA * SA* DB DZ+CA*CA* DT-sD RHO + SA*CA*l) B R D Z 
FS=C -BT>*C DRTDS>*GX*2. 0*PI *C 1 . 0/EMU) *RS**3 | 

EM U= 4 . 3 * P X * 1 . 3 E- 0 7 j 

FRHOA=FS*CA i 


FZACA=FS*SA 
FZ AP= FZACA 


FYAP=FRHOA*SIN( PHI A) ! 
FXAP=FRHOA*COSCPHI A) 'I 
FXA=F1*FXAP+F3*FVAP+F5*FZAP j 
FYA=<-F2>*FXAP+G1*FYAP+F6*FZAP 'j 
FZA=C-F4)*FYAP+F7*FZAP 

XC=F1*CX+G1*S)+F2*< y-F2*S> ] 


r'C = C -F3) *(X+G1*S>+G1*C Y-F2+S) + F4*< Z+G2*S> 


ZC = F5*(X+G1 *S>-F6*C y-F2*S)+F7*C Z+G2*S> 

zcc=zc ; 

RHOC=SQRTCXC*XC+rC*rC) j 

PHIC=ATAM2C VC/XC) j 

BPrll C= 0*0 ' ; 

CALL MAGC J/AM/EYC/RHOC/ZCC/BRHO/BZAC) j 

BRC 1 =BRHO 1 

IFCRHOC.LT. 1.0E-05) BRC 1=0.0 | 

BZC 1 =3ZAC j 

RHOC=RHOC+0. 0301 : 

CALL MAGC J/ AM/ ETC/ RHOC/ ZCC/BRHO/BZAC) : j 

BRC2=3RH0 \ 

BZC2=BZAC - ‘ 

D3DRH0= C BRC 2 -BRC 1 J/0.0001 % > 

RHOC=RHOC- 0 .0001 ’ 1 

ZCC=ZCC + 0. 3001 . . _ .. . ..... J 
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uoyyw 

01000 
010! 0 
01020 
01030 
01043 
01050 
01060 
01070 
01080 
01090 
01 100 
01 1 10 
01120 
01 130 
01 140 
01 1 50 
01160 
01 170 
01180 
01 190 
01200 
01210 
01220 
01230 
01240 
01250 
01260 
01270 
01280 
01290 
01300 
01310 
01320 
01330 
01340 
01350 
01360 
01370 
0138 0 
01390 
01400 
01410 
01420 
01430 
01440 
01453 
01460 
01470 
01480 
01490 
01500 
01510 
01520 
01530 
01 540 
01553 
01560 
01570 
01580 




CALL MAGC J, AM, EYC, RHOC, ZCC , BRHO, BZ AC ) 

BRC3=BRH0 
B7. C3 ---30 AC 

D3DZ = ( DZC3-BZC 1 )/0.0601 
DB HDZ- C B RC 3 - BRC 1 >/ 0 . 0 0 3 1 
DBDZR= C BZC2-BZC 1 ) /0. 000 1 
BT- 5QRTC BRC 1 +3RC 1 +BZC 1 -i-BZC 1 ) ! 

SA-BZC1/BT j 

CA-BRC 1 /3T ■ . j 

DBTDS=SA*CA*DBDZR+SA*SA*DBDZ+CA*CA*D3DRHO+SA*CA*DBRDZ 
FS-C -BT)*C DBTDS)*GX*2.0*PI*C 1 . 0/EMU) *RS**3 I 

FRHOC= FS*C A 

FZ ACC= FS*SA | 

FZCP=FZACC j 

FYCP=FRHOC*SIMC PHIC) | 

FX C P= FRf I OC *CO SC PHIC) . j 

FXC = F 1 * FXC P- F 3 * FYC P+ F 5 * FZ C P 
FYC = F 2 *FXCP-» G 1 *FYCP-F 6 +FZCP 
FZC = F 4 * F YC P+ F 7 * FZC P . 

XD=Y 

YD= C -F 7 ) *C X- F 7 *SQRTC 3 . 0 >*S>-F 4 *< Z+G 2 *S) 

ZD-( -F 4 )*C X-F 7 *SQRTC 3 . 0 ) *S) +F 7 *C Z+G 2 *S> 

ZDC--ZD 

RHOD- SORT C XD*XD+YD*YD) I 

PHID=ATAN2C YD,XD) •] 

BPH I D= 0 • 0 - ] 

CALL MAGC J, AM, E7D, RHOD, ZDC, BRHO, BZAC > ' j 

BHD 1 =8RH0 i 

IFCRHOD.LT. 1.0E-05) BRD1=0.0 I 

BZD! =BZAC • { 

RHOD=RROD-i . 0001 ’ 

CALL MAGC J, AM, EYD, RHOD, ZDC, BRHO, BZAC ) j 

BRD2=3RH0 j 

BZD2=BZAC j 

DBDRHO= C BRD2-BRD1 ) / • 000 1 i 

RHOD=RHOD-.0031 \ 

ZDC=ZDC+. 0001 ] 

CALL MAGC J, AM, EYD, RHOD, ZDC, BRHO, BZAC ) 1 

BRD3=BRH0 ] 

BZD3=BZAC 1 

D3DZ = C BZD3-BZD1 )/. 0001 ’ j 

DBRDZ=CBRD3-BRD1 )/0.0001 • j 

DBDZR=C BZD2-BZD1 )/0.0001 j 

BT-SQRTC BRD1 *BRD1 +BZB1 *BZD1 ) ] 

SA=BZDi/BT I 

CA=BRD1/BT ' j 

DBTDS= SA*C A+D3DZR+SA*SA+ DBDZ + CA*C A + DBDBHO + SA*C A*D8RPZ 
FS-C -BT)*C DBTDS) *GX*2« 0*PI *< 1 . 0/EMU) *RS**3 j 

FRHOD=FS*CA I 

FZAC D=FS*SA 1 

FZDP= FZACD j 

F7DP= FRHOD*S I X C PH ID) • ] 

FXDP=FRHOD*COSC PHI D) I 

FXD=C -F7 )*FYDP-F4*FZDP 
FYD= FXDP 

FZD= C - F 4 ) *FYDP+F 7 * FZDP 

X3-X j 

73= -Y 1 

Z 8 =-Z+ 0 . 25 *SQRTC 6 . 0 >*S j 
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r 01590 

RH0B~ SORT C XB*XD+ l'B*YB) 

' * ' 

i 01600 

PHIB=ATAN2C YB>XB) 

T 

01610 

‘ BPHIB=0.0 


[ 01620 

CALL MAGC J* AM* EYB* RHOB* ZBC * BRHO* BZAO 

\ 

; 01630 

BRB1-BBH0 

\ 

1 0.1640 

I FC RH OR <• L T • 1 • 0 !:> Q 5 ) BRB1=0.0 

1 

| 01650 

BZB 1 =BZ AG 

i 

\ 

i 01660 

RH0B= RHOB-i-0. 000 1 

s 

i 01670 

CALL MAGC J* AM> EYB* RHOB* ZBC* BRHO* BZAC > 


[. 01680 

BRB2=BRHO 


01690 

BZB2=BZAC 

l 

j 01700 

DBDRH0= ( BRB2-3RB 1 > /• 0001 


• 017 10 

RHOB=RHOB-.0001 

) 

j 01720 

ZBC=ZBC+» 0301 . • 

4 

1 

; 01730 

CALL MAGC J> AM* EY3* RHOB* ZBC* BRHO* BZAC ) 

i 

4 

! 01740 

BRB3=BRrlO 

. \ 

i 017 50 ' 

BZB3=BZAC 

, i 

1 01760 

DBDZ-CBZB3-BZB1 )/. 0001 


; 01770 

DBRDZ= C 3R33-BRB1 > /0. 000 1 

j 

i 01780 

DBDZ R- C BZB2-BZB 1 ) / 0 * 300 1 

] 

01790 

BT= SQRTC BFiB 1 *BR3 1 + BZ3 1 +P.ZD 1 ) 


01800 

SA=BZB 1 /3T 


018 13 

CA=BRB1/BT 

1 

! 01820 

DBTDS=SA*CA*DBDZR+ SA*SA*DBE»Z-»CA*CA*DBDRHO+S 

A*CA*DBRDZ 

! 018 3 0 

FS=< - 3 T ) * C DBTDS>*GX*2. 0*PI *< 1.0/EMU) *RS**3 

A 

1 01640 

F R H 0 B « V .S * C A 

1 

j 018 50 

FZ ACB-- FS * S A 

i 

t 

018 60 

FZBP= FZ AC 3 


I 01070 

FYBP«FRHOB*SirK PHIB> 

\ 

i 

1 01880 

FXBP=FRH08*COSC PHIB) 


; 01890 

FXB= FXBP 

i 

i 

01900 

F/B=-FYBP • 

1 

! 01910 

FZB=- FZBP 


01930 

FXC J ) " FXA+ FXB+ FXC Jr FX D 

i 

i 

[ 01940 

FTC J) ~ FYA+FYB+FTC+FYD 

j 

( 01950 . 

FZC J>= FZA+FZB+FZC+FZD 

1 

[ 01970 

59 CONTINUE 

\ 

t 

j 01980 

S UMX-0.0 

4 

i 

i 01990 

SUM Y= 0.0 

j 

02800 

SUMZ = 0 « 0 

i 

j 

1 02010 

DO 60 M=l*4 

\ 

! 02020 

SUMX= SUMX+FXC M ) 

. i 

i. 02033 

SUMY-SUMY+FYCM) 

i 

; 02040 

SUMZ-SUMZ+FZCM) 

i 

\ 

1 02050 

60 CONTINUE . 


‘ 02060 

SUMX=SUMX*9.0*1.0E+05 

•j 

i 02370 

S UM Y = S UM Y * 9 • 0 * 1 . 0E+05 

3 

1 

' 02083 

SUMZ = SUMZ*9. 0*1 • 0E+05 

< 

j. 02110 

FT0T= SQRTC SUMX+S UM X + S UM Y*S UM Y+SUMZ+S UM Z ) 

j 

! 02 1 50 

RETURN 

.j 

i 02160 

-END ■„ 



• 
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f 000 10 ' 

SUBROUTINE MAGC J* A* EYE* RHO, Z, BRHO, BZ ) 

■.I ... 

1 

00020 

DIMENSION EYEC 6> 

1 

■ - 1 

00030 

PI=3. 1415926 

) 

00040 

EMU=4.0*PI * 1 • 0E-37 

.j 

00050 

S 1=A*A >- RHO*RHO+Z*Z 

j 

00060 

S2= C A- RHO )*( A- RHO ) + Z*Z 

i 

i 

000 70 

S3-( A + RHO ) *( A+RHO) + Z*Z 

I 

1 00080 

S4=A*A-RH0*RH0-Z*Z 

i 

1 

i 03090 

EKS=4» 0*A*RHO/S3 

,’i 

| 00130 

EK=SGRTCEKS> 

j 

! ■ 00110 

Y=CELIC1,EK> 

< 

j 

, i 

00120 

TK=Y 

j- 00130 

Y=CELI ( 2, EK) 

'i 

1 

} 00140 

E=Y 

j 00150 

S5=-TK+< S1/S2)*E 

* 

I 00160 

S 6= Z / ( RH 0 * SORT (S3)) 

* 

; 00170 

BRHO=< ( EYEC J)*EMU>/C 2.0*PI ) )*S6*S5 

j 

! 00180 

S7= TK+C S4/S2 ) *E 


1 00190 

3Z=( ( EYEC J)*EMU) /C 2. 0*PI ) ) *C 1 • 0/SORTC S3) )*S7 


00200 

RETURN ' 

j 

00210 . 

END.- 

! 

i f 00010 

SUBROUTINE CNVRTFC EXP, FTP, FZP, FFXU, FFYU, FFZU, 

XMAGFF); 

! 00020 

A=SQRTC 1 . /3* ) 

« 

j 00030 

B= SQRTC 2 • /3 • ) 


! 00040 

FX= A * FX P +0.0 + FYP -B*FZ P 

1 

| 00050 

FY=0»0*FXP +1.0*FYP +0.0*FZP 

' 

00060 

FZ= B*FXP +0 • 0*FYP +A*FZP 


00070 

XMAGFF=SQRT( FX+*2+FY**2+FZ**2 > 


1 03080 

IF ( XMAGFF* EQ* 0* 0) GO TO 20 


| 00090 

* F UNIT VECTOR COMPONENTS 

i 

i 00100 

FFXU= FX/XMAGFF 


r. 00110 

FF YU= FY/XMAG FF 


j- 00120 

FF3U=FZ /XMAGFF 


i 00130 

GO TO 15 

i 

00140 

20 FFXU= 0 • 0 i FFYU= 0*0 ;FFZU=0.0 

j 

001 50 

15 CONTINUE 

j 

00160 

RE TURIN 

j 

... , 00170 

END. _ 

.. . ..... j 
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USFC - Form 3416*2 (February 1970) 


SLIST TETRA 


0 6/0 7/7?. 

0000 5 
00020 
00030 
00040 
000 50 
000 60 
000 70 

00125 

00126 
00 127 
00130 
00163 
00180 
002 10 

002 70 ’ 

00280 
00290 
00300 
00320 
00330 
00340 
00350 
00360 

003 70 
00380 
00390 
00400 
00405 
00410 
00420 
00430 
00440 
00450 
00 A 60 
00520 
00 530 
00540 
00545 
00550 
00590 
00 600 
00610 
00620 
00630 
00 640 
00650 
00660 
00670 


09*628 

DIMENSION EYAC4)* EYB C 4 > > EYC < 4 ) > EYD C 4 > > SK C 4 > > FX < 4 > > FY < 4 > > FZ ( 4 ) 
79 READsX^YjZ 

PRINT \ jXt Y tZ 

1 FORMAT (3X*2HX=F7.4#3X*2HY=F7.4*3X*2HZ=F7. 4) 

PRINT 2 

2 FORMAT (/> 

A=1 . 15 

EYA( J > = 1 « 0 ; EYA ( 2 ) =4 7 • 5;EYA(3)=0*8;EYA<4>=1 .0E-06 
EYRC 1 )=0 . 94; EYB(2 )=2. 41 EYB ( 3 ) =0 « 3 1 J EYBf 4) = 1 .0E-0 6 
EYC C 1 ) = 1 «82;EYC (2> = 1 • AS EYC < 3 > = 45.0; EYC (A )= 1 . 0E-06 
EYDC 1 ) = 4 5 . 0 ; EYD <2 > = 2 . 75; EYDC3>=2. 1 i;EYDC4)= 1 .0E-06 
SCM=3 .05 

READ; (SKCI>#I=1*4> 

RS = 0 . 5 

s=scm* .01 ;x=x*.0 1 ; Y--Y+.0 1 ; z=z*.0 1 

AM=A* . 0 1 
RS=RS*.01 
PI=3. 1415926 

F1=0.5*SQRTC3.0) .'*■ 

G1=C 1 .0/6.0 )*SGRTC3.0> 

F2=0 • 5 
F3=l. 0/6.0 

F4=C2.0/3.0)*SGRT(2.0) 

G2 = < 1 .0/12.0)*SGRTC6.0) 

F5=C 1 .0/3.0 )*SGRT (2.0) 

F6=( 1 .0/3.0 >*SGRT( 6.0) 

F 7=1. 0/3-0 
DO 59 J= 1 *4 

XA=F1*<X+G1*S>-F2*<Y+F2*S> 

YA=F3*<X+G1*S)+G1*(Y+F2*S>-F4*(Z+G2*S> 

ZA=F5*<X+G1*S)+F6*<Y+F2*S>+F7*(Z+G2*S> 

ZAC=ZA 

RHOA=SQRT <XA*XA+YA*YA) 

PHI A=AT AN2 C Y A > XA > 

BPHI A=0 . 0 

CALL MA G( AM* EYA*RHOA* ZAC* BRHO* BZAC ) 

BRA1 =BRHO 

I F (RHOA . LT . 1 «0E-0 5 ) BRA1=0.0 
BZA 1 =BZAC 

SIH=<EXP(2.0*SK(J> )=EXP(“SK( J) ) )/2.0 
SIH1 =<EXPCSK C J) )-EXPC-SK< J) ) >/2.0 

0= CS I H-S IN ( 2 .0*SK C J ) ) ) / C S I H 1 *S I H 1 +SI N < SK ( J ) ) *S I N(SK ( J > ) ) 

GX= 1 .0-(0. 75/SK(J) > *Q 
DRH0A= .000 1 
RH0A=RH0A+DRH0A 

CALL MAGC J, AM* EYA* RHOA* ZAC* BRHO « BZAC > 

BR A2=BRH0 

__BZA2=6ZAC _____ , 
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00680 
00 690 
no 70 ci 
00 7 1 0 
00700 
00 730 
00 735 
00 730 
0073 1 
00 742 
00 743 
00 744 
00745 
00 74 6 
00 74 7 
00 7 50 
00 760 
0 0 7 70 
00 780 
OO790 
00600 
006 10 
O0820 
00030 
OO3 70 
00 860 
00 890 
0 0 9 00 

or, ? i o 

00920 
00980 
OO990 
0 1 000 
0 100 5 
0 1 0 1 0 
0 1 0 50 
O 1 0 60 
0 1 OIO 
O 1080 
0 1 0 90 
0 1 1 00 
OHIO 
0 1 120 
0 1130 
0 1 1 40 
0 1 1 50 
0 1 151 
0 11 
0 1 1 53 
0 1 1 54 
0 1155 
0 1156 
0 1157 
...a.i_L Ad 


DBDrTiO=(BRA2-BRA1 J/DRHOA 

RHOA=RHOA-DRHOA / 

Z AC = Z AC ♦ . 000 1 

GALL MA G( J* AM* F.Y A* RHOA* ZAC »BRHO*BZAC> 

BRA3 =8RH0 
BZA3=BZAC 

F.MU = 4 .0*P I * 1 • 0E“0 7 
D 0 D Z = ( B Z A 3 - D Z A 1 >/0.0001 
DBRDZ=(BRA3-BRA| )/0.0001 
DBDZR= U3ZA2-8ZA1 )/0.000l 
BT = SQRT (BR A 1 *BR A 1 +BZA 1 *BZA 1 ) 

SA-nZA 1 /BT 
CA=BRA 1 /BT 

DBTDS = S A + C A* DBDZFC + SA + S A+DRDZ+CA + CA*L)BDRHO + S A*C A + DBRDZ 

r 51 = (-RT > + CDBTDS )*GX*2 .0*PI* ( 1 .0/EMU>*RS**3 

E;mii = 4.0*PI* I « 0 F. - 0 7 

FPHOA=FS*CA 

FZACA=FS*SA 

FZAP--FZACA 

F YAP = FR HO A *3 I N ( P H I A ) 

FXAP = FRHOA*COS (PHI A) 

F X A = F l * FXAP+F 3*F YAP+F 5*FZ AP 
FYA-- ( -F2 ) + FXAP+ G 1 *F YAP + F6*FZAP 
F ZA= C-F4 >*KYAP+F 7*FZAP 
X.C=F 1*(X+G1 +S >+F2*(Y-F2*S> 

YC=(-F3 >* (X+GMS ) + Gl * ( Y- F2*5 ) + F4* ( Z+ G2*S > 

ZG-F5* (X+01 *S)-F6*(Y-F2*S > + F7* (Z+G2*S> 
zoc-zc 

Fv HOC “3 CRT (XC*XC + YC*YC ) 

PI I C = AT AN 2 ( YC * XC ) 

BPHIC=0.0 

CALL MAC ( J* AM* EYC*RHOC* ZCC* BRHO, BZAC > 

PRC I =RRHO 

IFCRHOC.LT. 1 • 0 E - 0 5 ) BRC 1=0.0 

H?:0 1 sc ZAC 

RROC = .RHOC+0 • 000 1 

CALL MAG< J* AM, E YC * RHOC * ZCC * BRHO * BZ AC ) 

His C2 = BR HO 
L<ZC2 = BZAC 

DP l)R HO = (BRC2-BRC 1 )/0 .000 1 
W HOC =R. HOC -O.000 1 
Z L L = Z o C + 0 . O' 0 vj 1 

CALL MAGC J* AM* E YC * RHOC * 7.CC * BRHO* BZ AC ) 

fH-;C3 = ERH0 

PZC3=BZAC 

Dr')Z= (BZC3-BZC 1 )/0 .000 1 
D ; RDZ=(PRC3-BRC 1 )/0 .090 1 
PPDZR- CPZC2 -HZC 1 )/0 .990 1 
BT--SGRT (BRC 1 *BRC 1 ♦BZC 1 *ezC 1 ) 

S A = 5 Z C l /BT 
CA=BRC 1 /BT 

DBT0S=SA*CA*DBDZR+SA*SA*DBDZ+CA*CA*DBDRH0+SA+CA*D3RDZ 

FS= C -BT >* (DOTDS ) *GX + 2 .9*PI * ( 1 . 0/ EMU ) *RS ** 3 

CBunr«EC*<«A. . — — 






01160 FRHOC=FS*CA 

01170 FZACC=FS*SA 

01180 FZCP=FZACC 

01190 FYCP=FRHOC*SIN(PHIC) 

01200 FXCP=FRHOC*COS (PHI C > 

01210 FXC=F1*FXCP-F3*FYCP+F5*FZCP 

01220 FYC=F2*FXCP+G1*FYCP-F6*FZCP 

01230 FZC=F4*FYCP+F7*FZCP 

01270 XD=Y 

01280 YD=C-F7)*CX-F7*SQRTC3.0>*S>-F4*(Z+G2*S> 

01290 ZD=C-F4)*CX-F7*SGRT<3.0)*S)+F7*(Z+G2*S> 

01300 ZDC =ZD 

01310 RHOD=SQRT C XD*XD+ YD* YD ) 

01320 x PHID=ATAN2 CYDjXD) 

01380 BPHID=3.0 

01390 _ CALL MAG( J* AM* EYD.*RHOD.» ZDC * BRHO* BZAC > 

01400 BRD 1 =BRHO 

01405 IFCRHODoLT. 1 .0E-05) BRD1=0.0 / 

01410 BZD 1 =BZAC 

01450 RHOD=RHOD+.000 1 

01460 CALL MAGC J, AM* EYD* RHOD> ZDC j BRHO* BZAC > 

01470 BRD2=BRH0 

01480 BZD2=BZAC 

01490 DBDRHO= CBRD2-BRD1 )7 • 000 1 

01500 RHOD=RHOD-.000 1 

01510 ZDC=ZDC+.0001 

01520 CALL MAGC J * AM* EYD* RHOD* ZDC * 8RH0* BZAC > 

01530 BRD3=BRH0 

01540 8ZD3=8ZAC 

01550 DBDZ=(BZD3-BZD1 J/.0001 

01551 DBRDZ= ( BRD3 -BRD 1 )/0*0001 

01552 DBDZR = ( BZD2-BZD 1 J/0.0001 

01553 BT=SQRT CBRD 1 *BRD 1+BZD1*BZD1 > . 

01554 SA=BZD1/BT 

01555 CA=BRD1/BT 

01556 . DBTDS = SA*CA*DBDZR + SA*SA*DBDZ-i-CA+CA*DBDRHO + SA*CA*DBRDZ 

01557 FS= ( -BT )*(D3TDS)*GX*2»0*PJ*( 1 . 0/EMU > *RS**3 

01560 FRHOD=FS*CA 

01570 FZACD=FS+SA 

01580 FZDP=FZACD 

01590 FYDP=FRHOD*SJN<PHID> 

01600 FXDP=FRHOD*COS(PHID> 

01610 FXD=C-F7)*FYDP-F4*FZDP 

01620 FYD=FXDP 

01630 FZD=C-F4)*FYDP+F7*FZDP 

01670 XB=X 

01680 YB= - Y 

01690 ZB=-Z+0.25*SGRT< 6.0)*S 

01700 ZBC=ZB 

01710 RHOB=SGRT CXB*X8+YB*YB ) 

01720 PHI B = AT AN2 ( YB* XB > 

01780 EPHIB=0.0 

01790 CALL MAGC J* AM* EYB,RHOB* ZEC* BRHO* BZAC > 

-Cu t oaa odo « »pp uq . . . „ , , __ - 
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0 1 800 


BRB 1 =BRHO 

01805 


J F (RHOB « LT • 1 • 0E-0 5 > BRB1=0.0 

01810 


BZB 1 = BZAC 

0 1 860 


Rh'OB = RHOB+0 .000 J 

01870 


CALL MAG< J* AM* EYB# RKOB> ZBC » BRHO* BZAC ) 

• 01880 


BRB2=BRH0 

0 1 890 


BZB2=BZAC 

01900 


DBDRH0=(BR82-BRB 1 >/.000 1 

01910 


RHOB =R HOB- .000 1 

01920 


ZBC=ZBC+ .000 1 

01930 


CALL MAG<J.»AM.>EYB*RHOB*ZBC*BRHQ*BZAC> 

01940 


BRB3=BRH0 

0 1950 


BZB3=BZAC 

01960 . 


DBDZ= f BZB3-BZB 1 >/.000 1 

01961 


DBRDZ=<BRB3-BRB1 )/0.0001 

01962 


DBDZR = (BZB2-BZB 1 )/0 . 000 1 

01963 


BT=SQRTCBRB1*BRB 1+BZB1*BZBI > 

0 1964 


SA=BZB1/BT 

01965 


CA=BRB1/BT 

0 1966 


DBTDS=SA*CA*DBDZR+SA*SA*DBDZ+CA*CA*DBDRHO+SA*CA*DBRDZ 

01967 


FS=C-BT)*CDBTDS)*GX*2.0*PI*C1.0/EMU)*RS**3 

01970 


FRHOB=FS*CA 

0 1 980 


FZACB=FS*SA 

0 1 990 


FZBP=FZACB 

02000 


FYBP=FRHOB*SINfPHIB> 

02010 


FXBP=FRMOB+COS (PHIB ) 

02020 


F XB=FXBP 

02030 


FYE=-FYBP 

920 40 


FZB=-FZBP 

020 60 


PRINT 2 

02080 


FXf J)=FXA+FXB+FXC+FXD 

02090 


FY<J)=FYA+FYB+FYC+FYD 

02 100 


FZ< J)=FZA+FZ8+FZC+FZD 

02 120 

50 

FORMAT <3 X> 3HFX=E1 4.6*2X*3HFY=£14.6>2Xj3HFZ=E14.6> 

02 125 

59 

CONTINUE 

02 1,30 


SUMX-0 . 0 

02 1 40 


SUMY =0.0 

02 150 


SUMZ=0 . 0 

02 1 60 


DO 60 M = 1 r A 

02 1 70 


SUMX = SUMX + FX CM > 

02 1 80 


SUMY=SUMY+F YCM) 

02190 


SUMZ=SUMZ+FZCM) 

02200 

60 

CONTINUE 

02210 


SUMX = Sl'MX*9.0*l *0E+0'5 

02220 


SUMY =SUMY*9 • 0* 1 . 0E+05 

02230 


SU!<Z = Sl!MZ*9.0* 1 .0E+0 5 

02240 


PRINT : "COMPUTE TOTAL FORCE COMPONENTS" 

022 50 


PR I NT 50 * S UM X t S !JM Y > S UM Z 

02260 


F l OT = SGRT C SUMX*SUUX+3UMY*SUMY+SUM7*SUMZ ) 

02270 


PRINT 70> FTOT 

02280 

70 

F0RMATC3X, 12HT0TAL F0RCE=E14.6, 1X,5HDYNES> 

. 02282 


PR I NT : " GO ,TO NEXT CASE" 

02283 

<aoo a* 


PRINT 2 , . 
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02284 

02290 


GO TO 79 
END 


READY 


1ST CUBE 
06/07/72 09 « 32 5 




00010 
00020 & 
00030 
00040 
000 50 
000 60 
00070 
00080 

00090 

00091 
00100. 
00101 
00110 
001 11 
00120 
00121 

00130 

00131 
00 1 40 
00141 

00150 

00151 
00 1 60 
00 1 70 
00180 
00 1 90 
00200 


EYC C 1 ) = 1 
EYC ( 5 ) = 1 
EYD < 1 > = 1 


D , MENS ION EYA<6>.EYB<6>.EYC<6>.EYD<6>,EYE<6>,EYF<6>, 

SKC6)j FX<6>*FY<6)*FZ<6> 

79 READ:X*Y*Z 

, FORMAT ( 3X»2HX=F7« A* 3X»2HY = F 7 •4#3X>2HZ S F7»4) 

PRINT 2 
2 FORMAT </> 

EYA( !)-. •OE-06 J EYA<2 )a1»0E*06JEYAC3) s A5»0*EYA(4>-l*0E - 06 

EYB M ) = 1 *0E-3 6 J EYB <2 > * 1 *GE~06JEYB<3) = 1 • 0E*06J EYB ( A ) = I «0E“0< 
E YE ( 5)M.0E:06;EYB<6>=1.®E-06 j eyc (3>s1 . 0 e-06;EYC(4.= I.0E-0« 

tIlJV j , .0E-06JEYD(2)=A5^0JEYD(3>= 1«0E*06J £YD<A) s 1 • 0E-06 

SSn!:!:SE-S«wr«U!:SE-S6,EYE«3».|.K-B6.EYEC«.l.»M( 

^r«r.Ui”;wF«»-? > -« : «^F<3»...eE-«6.EYF<4,...«-.6 
SKC5)=1 .0E-06*SKC6>=l*0E-06 

SCM=3 • 33 

S=SC0*0.01JX=X*0.01iY=Y*0*0 1 iZ=Z*0' 01 

AM=A*0«01 
RS = RS*0 *0 1 

_,t»i «.o. — i. A I.Q.Q9..A — — »- — -- — ^ 
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00210 

00220 

00230 

00240 

00250 

00260 

002 70 - 

00280 
00290 

003 00 
00310 
00320 
00330 
00340 
00350 
00360 
00370 
00380 
00390 
00400 
00410 
00420 
00430 
00440 
00450 
604 60 

004 70 
00480 
03490 
00500 
00510 
00 520 
00530 
00540 
00550 
00 560 
00570 
00 580 
00 590 
00 600 
00610 
00620 
00630 
00640 
00 650 
00 660 
00670 
0.0 680 
00 690 
00 700 
00710 
00 720 
00 730 


PI =3 • 1 4 1 592 6 
no 59 J= 1 * 6 
XA = Z 
YA = Y 

ZA=S/2.0-X 
ZAC = ZA 

RHOA = SQRT < XA* XA+ YA* YA > 

'PHI A=AT AN2 ( YA# XA) 

BPHI A=0 .0 

CALL MA Gf J# AM# £YA»RHOA» ZAC# BRHQ # BZAC ) 

BR A 1 =BRHO 

IFfRHOA.LT. 1 .0E-05) BRA1=0.0 
BZA 1 =BZ AC 

S I H= ( EXPC2 .0*SK< J ) ) ~ EXP C -SK < J ) ) )/2 .0 
SI HI =<EXP(SK(J> )-EXP<-SK<J> ) )/2.0 

Q= <S I H-SIN<2 «0*SK( J > ) >/ <SI HI *SIH 1 +SIN (SK< J) >*SINCSK ( J> > > 
GX=1 .0-C0. 75/SKCJ) )*0 
RHOA=RHOA+0 .000 1 

CALL MAGf Jj AM* EYA#RHOA* ZAC # BRHO# BZAC ) 

BR A2 =BR HO 
BZA2=BZAC 

DBDRHO= (BRA2-BRA1 )/0 .000 1 
R HO A =R HO A -0 .000 1 
ZAC = ZAC +0.000 1 

CALL MAG< J» AM,EYA.,RHOA#ZAC#BRHO#BZAC> 

BR A3=BRH0 
BZA3=BZAC 

EMI_I=4 . 0 *P I * 1 .0E-07 
Dr>DZ = C 3ZA3-BZA 1 >/0 *009 1 
DBRDZ= (BRA3-SR A 1 )/0.0001 
DBDZR= (BZA2-BZA 1 J/0.000 l 
BT=SQRT<BRA1*BRA1+BZA1*BZA1 ) 

SA=BZA1/BT 
C A=5R A 1 /BT 

DBTDS=SA*CA*DBDZR+SA*SA*DBDZ+CA*CA*DBDRH0+SA*CA*D6RDZ 

FS = ( -BT ) * CDBTDS >*GX*2»0*PI*< 1 • 0/EMU > *RS**3 

FR HOA=FS*C A 

FZACA=FS*SA 

FZAP=FZACA 

FY AP=FRHOA*S I N < PHI A ) 

FXAP=FRHOA*COS (PH I A ) 

FXA=-FZAP- 
FYA=FYAP 
FZA=FXAP 
XB = X 
YB--Z 

ZB = Y +S/2 .0 

7nr. = 7B 

RHOB=SGRT CXB*X5+YB*YB ) 

PHI B = AT AN2 ( YB» XB ) 

BPHI B=0 • 0 

CALL MAG< J# AM# EYB #RHOB# ZBC # BRHO » BZAC ) 

BRB 1 =BRhO ... _ — — * 



00T40 
00750 
00760 
00770 
00780 
00790 
00800 
00810 
00820 
00830 
00840 
00850 
00860 
00870 
00880 
00890 
00900 
00910 
00920 
00930 
00940 

009 50 
00960 
00970 
00980 
00990 
01000 
01010 
01020 

' 01030 
01040 
01050 

010 60 
0 10 70 
0 1080 
0 1090 
01100 
0 1110 
0 1 120 
0 1130 
01140 
01150 
0 1 1 60 
01170 
0 1180 
0 1190 
01200 
012 10 
0 1220 
0 1230 
01240 
0 12 50 
01260 

L 0JLQ 7CI _ 


" IF CRHOB.LT. 1 .0E-05) BRB1=0.0 ' l ; " ■ " 

BZB1 =BZAC 
RHOB=RHOB+0 .000 1 

CALL MAGC J* AM* EYB*RHOB* ZBC* BRHO*BZAC > 

BRB2=BRH0 

BZB2=BZAC 

DBDRHO= C8RB2-BRB 1 )/0.0001 
RHOB=RHOB-0 .000 1 
ZBC = ZBC +0 *0001 

CALL MAGCJ*AM*EYB*RHOB*ZBC* BRHO* BZAC > 

BRB3=BRHO 

BZB3=BZAC 

DBDZ=(BZB3-BZB1 )/0.0001 
DBRDZ= CBRB3-BRB 1 J/0.000! 

DBDZR= CBZB2 ~BZB 1 J/0.0001 
BT =SGRT CBRB 1 *BRB 1 +BZB 1*BZB 1 ) 

SA=BZB1/BT 

CA=BRB1/8T 

DBTDS=SA*CA*DBDZR+SA*SA*DBDZ+CA*CA*DBDRHO+SA*CA*DBRDZ 

FS=C-RT)*CDSTDS)*GX*2.0*PI*C 1 . 0/ EMU ) *RS**3 

FRHOB=FS*CA 

FZAC8=FS*SA 

FZBP=FZACB 

FYBP=FRHOB*SI NCPHIB) 

FX8P = FRHOB*COS C PHI B ) 

FXB=FXBP 
FYB=FZBP 
FZB=-FYBP 
XC = -Z 
YC=Y 

ZC=X+S/2.0 

ZCC=ZC 

RHOC=SQRT CXC*XC+YC*YC ) 

PHI C =AT AN2 C YC * XC ) 

BPH IC=0 .0 

CALL MAGC J* AM* EYC*RHOC* ZCC/BRHO* BZAC ) 

BRC 1 =BRHO 

IFCRHOC.LT. 1 .0E-05) BRC 1 =0 . 0 
BZC 1 =BZAC 
RHOC=RHOC+0.000 1 

CALL MAGC J, AM* EYC*RHOC* ZCC* BRHO* BZAC) 

BRC2=BRH0 

BZC2=BZAC 

DBDRHO= CBRC2-BRC 1 )/0.0001 
RHOC=RHOC-0 .000 1 
ZCC=ZCC+0 .000 1 

CALL MAGC J* AM* EYC*R HOC* ZCC* BRHO* BZAC) 

BRC3=BRH0 

BZC3=BZAC 

DBDZ=CBZC3-BZC1 )/0.0001 
D8RDZ=CBRC3“BRC1 )/0.0001 
DBDZR= (BZC2-BZC 1 )/0.0001 . 

BT = SQRTCBRC 1 *BRC 1 +BZC 1 *BZC 1 > 

1 /QT . 
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012 70 
01280 
01290 
01300 
01310 
01320 
01330 
01340 
0 1 3 50 
01360 
01370 
01380 
01390 
01400 
01410 
01420 
01430 
01440 
01450 

0 1460 
01470 
01480 
01490 

01 500 
01510 

0 1 520 
01530 
01540 
01550 

01 560 
01 570 
01580 
01590 
01600 
01610 
01620 
0 1 630 
0 1 640 
0 1 650 
0 1 660 
0 1670 
0 1 680 
0 1 690 
0 1 700 
01710 
0 1 720 
0 1 730 
01740 

0 1 750 

01 760 
0 1 770 
0 1 780 
0 1 790 


SA=BZC1/BT 

CA=BRC1/BT 

dbtds=sa*ca*dbdzr+sa*sa*dbdz+ca*ca*dbdrho+sa*ca*dbrdz 

FS=(-BT >*(DBTDS>*GX*2.0*PI *( 1 .0/EMU>*RS**3 

FRHOC=FS*CA 

FZACC=FS*SA 

FZCP=FZACC 

FYCP=FRHOC*SI NCPHIC ) 

FXCP=FRHOC*COS (PHIC > 

FXC = FZCP 
FYC=FYCP 
FZC=-FXCP 
XD = X 
YD = Z 

ZD=S/2.0-Y. 

ZDC=ZD 

RhOD=SORT(XD*XD+YD*YD> 

PHID=ATAN2(YD*XD> 

BPH I D=0 • 0 

CALL MAG(J,AM*EYD>RHOD>ZDC>BRHO>BZAC) 

BRD 1 =8RM0 • 

IFtRHOD.LT. 1 .0E-05) BRD1=0.0 
BZD 1 =BZAC 
RHOD=RHOD+0.000 1 

CALL MAG ( Jj AM* EYD* RHOD* ZDC* BRHO* BZAC ) 

BRD2=BRH0 
BZD2 =BZAC 

DBDRHO- (BRD2-BRD 1 >/0.000 1 
RHOD = RHOD-0 «000 1 
ZDC=ZDC+0.0001 

CALL MA Gf J# AMj EYD* RHOD* ZDC * BRHO* BZAC ) 

BRD3=BRH0 
BZD3 =BZAC 

DBDZ=(BZD3-BZD1 )/0.0001 
DBRDZ= (BRD3-BRD 1 )/0.0001 
DBDZR=(BZD2-BZD1 )/0.000 1 
BT=SQRT (BRD 1 *BRD 1 +8ZD 1 *BZD 1 ) 

SA=BZD 1 /BT 
CA=BRD1/BT 

BBTDS = SA*CA*DBDZR + SA*SA*DBDZ+CA*CA*DBDRh'0 + SA*CA*DBRDZ 

FS=C-BT)*(DBTDS)*GX*2.0*PI*( 1 . 0/EMU >*RS**3 

FRHOD=FS*CA 

FZACD=FS*SA 

FZDP=FZACD 

FYDP=FRHOD*SIN(PHID) 

FXDP=FRHOD*COS(PHID> 

FXD= FXDP 
FYD=-FZDP 
FZD=FYDP 
XE = X 
YE = -Y 

ZE=S/2.0-Z 

ZEC=ZE _ 
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1 


0 I 800 
01810 
01820 
0 1830 
01840 
0 1 850 
0 1 860 
01870 
01880 
01890 
01900 
01910 
01920 
1 01930 

01940 
01950 
01960 
01970 
01980 
01990 
02000 
02010 
02020 
02030 
02040 
02050 
02060 
02070 
02080 
02090 
02 100 
02 110 
02 120 
02 1 30 
02 140 
02 1 50 
02 1 60 
02 1 70 
02180 
02 1 90 
02200 
022 10 
. 02220 
02230 
02240 
022 50 
02260 

022 70 
022 80 
02290 
02300 

023 10 
02320 


'RHOE=S0RT<XE*XE+YE*YE> i " 

PHI E=ATAN2 C YE* XE> 

BPHIE=0.0 

CALL MAGCJ*AM*EYE*RH0E*ZEC*BRH0>BZAC5 
BRE1 =BRHO 

JFCRHOE.LT. 1 .0E-05) BRE1=0.0 

BZE 1 =BZAC 

RHOE=RHOE+0.0001 

CALL MAG( Jj AM*EYE*RHOE* ZEC* BRHO*BZAC> 

BRE2=BRH0 
BZE2 =BZAC 

DBDRH0=CBRE2-BRE1 5/0.000 1 
RHOE=RHOE-0 .000 1 
ZEC=ZEC+0 .000 1 

CALL MAG C J # AM# EYE*RHOE*ZEC*BRHO* BZAC ) 

BF}E3=BRHO 

BZE3=BZAC 

DBDZ=CBZE3-BZE1 5/0 .000 1 
DBRDZ=CBRE3-BR£1 5/0.0001 
DBDZR=C6ZE2-BZE1 5/0.0001 
6T = SQRT CBRE 1 *BR£ 1+BZEJ+BZE1 ) 

SA=BZE1/BT 

CA=BRE1/BT * 

DBTDS=SA*C A+DBDZR +SA*S A*DBDZ+C A*CA*DBDRHO+S A*CA*DBRDZ 

FS=C-BT5*CDBTDS5*GX*2.0*PI*C 1 . 0/ EMU 5 *RS**3 

FRHOE=FS*CA 

FZACE=FS*SA 

FZEP=FZACE 

FYEP=FkHOE*SINCPHIE5 ' • 

FXEP=FRHOE*COS CPHI E5 

FXE=FXEP ' 

FYE=-FYEP 
FZE= -FZEP 
XF = X 
YF=Y 

ZF=Z+S/2.0 

ZFC=ZF 

RHOF=SGRTCXF*XF+YF*YF 5 
PHIF=ATAN2 CYF* XF5 
BPHI F=0 .0 

CALL MAGC J t AM* EYF * RHO F* ZFC * BRHO* BZAC ) 

BRF 1 =BRHO 

IFCRHOF.LT. 1 .0E-05) BRF 1 =0 . 0 
BZF 1 =BZAC 
RHOF=RHOF + 0 .000 1 

CALL MAGC J, AM* EYF* RHOF, ZFC* BRHO* BZAC 5 

BRF2=BRH0 

BZ F2 = BZAC 

DBDRH0=CBRF2-BRF1 5/0.000 1 
RHOF=RHOF-0.000 1 
ZFC = ZFC +0.000 1 

CALL MAGC J,AM*EYF*RHOF* ZFC* BRHO* BZAC 5 
BRF3=BRH0 

R m = R 7 Aft. . : - . . . ... 
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02330 1 ’ 


"STF3 — ;u,J " - ‘"-r-r-. ~ 

02340 


DBDZ=<BZF3-BZF1 5/0.000 1 

02350 


DBRDZ=(BRF3-BRF1 5/0.0001 

02360 


D6DZR= CBZF2-BZF 1 5/0.0001 

02370 


BT=SQRT<BRF1*BRF1+BZF1*BZF15 

02380 


SA=BZF1/BT 

02390 


CA=BRF1/BT 

02400 


DBTDS=SA*CA*DBDZR+SA*SA*DBDZ+CA*CA*DBDRHO+SA*CA*DBRDZ 

024 10 


FS = ( -BT 5 =MDBTDS 5 * GX*2 • 0*P I * f 1 • 0/ EMU 5 *RS**3 

02420 


FRHOF=FS*CA 

02430 


FZACF =FS*SA 

02440 


FZFP = FZ AC F - 

02450 


FYFP=FRH0F*SIN<PHIF5 

02460 


FXFP=FRH0F*C0S(PHIF5 

024 70 


FXF = FXFP 

02480 


FYF = FYFP 

02490 


FZF = FZFP 

02 500 


PRINT 2 

02510 


FXf J5=FXA+FXB+FXC+FXD+FXE+FXF 

02 520 


FYC J5=FYA+FYB+FYC+FYD+FYE+FYF 

02530 


FZ<J5=FZA+FZB+FZC+FZD+FZE+FZF 

02540 

50 

FORMAT <3X* 3HFX = E1 4. 6 > 2X* 3HFY=E14.6>2X>3HFZ = E14* 65 

02 550 

59 

CONTINUE 

02560 


SUMX=0 .0 

02570 


SUMY =0.0 

02 580 


SUMZ=0 • 0 

02590 


DO 60 M= 1*6 

02 600 


SUMX=SUMX+FXCM5 

02610 


SUMY=SUMY+FYCM5 

02 620 


SlJMZ = SUMZ + FZfM5 

02630 

60 

CONTINUE 

02 640 


SUMX=SUMX*9.0* 1 .0E+05 

02650 


SUMY=SUMY*9.0* 1 .0E+05 

02660 


SUMZ=SUMZ*9.0* 1 .0E+05 

02670 


PR I NT : "COMPUTE TOTAL FORCE COMPONENTS” 

02 680 


PRINT 50> SUMX* SUMY > SUMZ 

02690 


FTOT = SQRT< SUMX*SUMX+SUMY*SUMY+SUMZ*SUMZ > 

02 700 


PRINT 70 » FTOT 

02 710 

70 

FORMAT ( 3X j 12HT0TAL F0RCE=E14.6* 1X>5HDYNES5 

02 720 


PR I NT : ” GO TO NEXT CASE” 

02730 


PRINT 2 

02 740 


GO TO 79 i / 

02750 


END 


READY 
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APPENDIX C 


UNITS AND CONVERSIONS 

For magnetic field computations, CGS and gaussian units are most con- 
venient,, In this system of units based on tne centimeter, gram and second 
the magnetic field strength H in vacuum (or air) along the axis of a single 

2t 7 a 

turn coil of radius a is given by H = — — — i , where H is the field strength in 

r 

oersteds and r is the slant range to the coil winding, a and r are both 
measured in centimeters, and i is the current flowing measured in ab- 
amperes (10 amperes = 1 abampere). All of the coil configurations con- 
sidered in this study can be resolved into a summation of such single turn 
elements. In a vacuum the field strength H and the magnetic induction have, 

in CGS units, the same numerical value, i.e. , B = H. B is measured in gaus 

— * 1 -7 -* 

The force on a conductor carrying a current density j is given by — j x H 
dynes per cubic centimeter, where j is measured in abamperes per square 
centimer and c is the speed of light measured in centimeters per second. 

In MKS units, which are most convenient in discussing practical 

electrical circuitry and which are also becoming more widely used in all 

phases of electrical theory, the field strength H due to a single current 

9 . 

carrying coil turn of radius a is given by H = — — i along the axis of the 

2r 

winding at a slant distance r. The coil radius a and the slant range r must 
be expressed in meters and the current i in amperes. The unit of field 



strength H is the Maxwell (or ampere turn per meter). In this set of units 
B and H have a different numerical value in a vacuum. In MKS units the 
vacuum relation between B and H is given by B = 477 • 10 H. B is measured 
in webers per square meter. 


The conversion between these units is given in the brief list which 
should suffice for the limited scope of the present report. 


CGS-Gauss 


1 gauss 


1 oersted = 
1 abampere 
1 cm 



MKS 

in' 4 ' 2 

10 w-m 

Maxwells =79.5 Maxwells 

= 10 amperes 

10" 2 m 


C-2 



